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INTRODUCTION 

The interaction between the solar wind and the inter- 
stellar medium has been a topic of interest for many 
years [Davis, 1955, 1962; Parker, 1961, 1963; Axford 
etal, 1963] , but progress has been rather slow largely 
because of the paucity of relevant observations. 
Recently, however, the situation has changed as a result 
of the clear evidence for the penetration of interstellar 
gas into the inner solar system obtained from OGO 5 
[ Thomas and Krassa, 1971; Bertaux and Blamont, 
1971]. Furthermore, the prospect that there will be 
space probes traveling toward the outer regions of the 
solar system within a few years (see fig. 1) has created a 
need for a better understanding of the subject so that 
the most effective use can be made of the missions. 

This review is a survey of the work, both published 
and unpublished, that has been carried out up to the 
time of writing. To the best of our knowledge, the 
bibliography is complete up to mid-1971. First, the 
expected characteristics of the solar wind, extrapolated 
from the vicinity of the earth are described, and several 
simple models are examined for the interaction of the 
solar wind with the interstellar plasma and magnetic 
field. The possibility that there is a substantial neutral 
component of the interstellar gas is ignored. In later sec- 
tions, we consider various aspects of the penetration of 
neutral interstellar gas into the solar wind, and describe 
the dynamic effects of the neutral gas on the solar wind 
(to the extent that they are understood). Finally, we 
discuss problems associated with the interaction of 
cosmic rays with the solar wind. 

We briefly summarize here the present status of 
knowledge concerning the properties of the interstellar 
medium in the vicinity of the solar system [see, e.g., 
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Figure 1 Spacecraft trajectories for the proposed 
“Grand Tour” missions to Jupiter-Satum-Pluto 
(1976-1977), and Jupiter- Uranus-Nep tune (1979). The 
trajectories carry the spacecraft in the general direction 
of the projection of the solar apex into the ecliptic 
plane, the direction in which the distance to the solar 
wind termination is a minimum . The final portion of the 
trajectories (dashed) can he varied, since it depends on 
the last planetary encounter . 

Allen, 1963] . The sun is located on the inner edge of the 
local (Orion) spiral arm at a distance of ~8 parsec north 
of the plane of the galaxy. The solar motion relative to 
neighboring stars is 20 km sec" 1 in the direction 
a ~ 271°, 5 .** +30°, or l n ~ST,b n ** +22° (the “solar 
apex”). The nearest star is more than 1 parsec away, and 
the mean density of stars in the solar neighborhood is 
~0.057 solar masses per cubic parsec. Since the region 
containing the solar wind (the “heliosphere”) has dimen- 
sions that are at most of the order of 10 2 -10 3 AU 
(<10“ 2 parsec), and it seems reasonable to assume that 
stellar winds associated with nearby stars are not signifi- 
cantly more intense than the solar wind, we can be fairly 
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sure that the solar wind interacts directly with the inter- 
stellar medium and that there is no overlapping of stellar 
winds in our vicinity. It should be noted that the sun 
appears to be on the fringes of the Gum nebula; how- 
ever, according to Brandt et al [1971] it lies some 
60 parsecs outside and the local interstellar medium may 
not be directly affected by the presence of the nebula. 

The mean interstellar electron density n e in the disc 
of the galaxy can be determined in two ways: (1) by 
measurements of low frequency radio wave absorption 
[EUis and Hamilton , 1966; Smith , 1965; Bridle , 1969; 
Bridle and Venugopal, 1969; Alexander et al, 1969] ; 
and (2) by measurements of the frequency dispersion of 
pulsar emissions [Pilkington et al, 1968; Lyne and 
Rickett, \96%\ Habing and Pottasch, 1968 \Davies, 1969; 
Prentice and Ter Haar, 1969 a,b; Davidson and Terzian, 
1969; Mills, 1969; Gould, \91l\Grewingand Walmsley, 
1971]. The first method yields the quantity (rc e 2 7y~ 3/2 >/, 
where l is the distance in the galaxy along the line of 
sight and T e is the electron temperature; this tends to 
overemphasize the contributions from denser regions 
[Gould, 1971]. The second method yields <n e )l and thus 
permit, us to estimate <n e ) directly if we know the dis- 
tance to the pulsar by independent means (as in the case 
of the Vela and Crab pulsars), or if an estimate can be 
made of the neutral hydrogen content in the line of sight 
(e.g., from 21-cm absorption of the pulsar emission, as in 
the case of CP0329) so that a comparison can be made 
with 21-cm emission in the same direction. In fact, the 
interpretation of the observations is not easy because of 
the confusion introduced by dense ¥BI regions, and the 
structure associated with the spiral arms of the galaxy. 
According to Davidson and Terzian [1969] the mean 
electron density in the disc of the galaxy is 
0.03-0.10 cm" 3 . However, the result is model dependent, 
and according to Gould [1971] , who notes that the path 
to the Crab pulsar crosses the interarm region where the 
density could be low, the local mean electron density in 
the galactic plane is — 0.12 cm" 3 . 

Evidence for the presence of a large-scale interstellar 
magnetic field arises from observations of (1) the polari- 
zation of starlight [Davis and Greenstein , 1951; 

Mathewson, 1968] ; (2) the polarization of galactic non- 
thermal radio emission [Mathewson and Milne, 1965]; 
(3) the Zeeman effect at 21-cm wavelength [Verschuur, 
1968, 1969 a,b,c, 1970; Davies etal, 1968]; (4) the 
Faraday rotation of polarized extragalactic radio sources 
[Gardner and Whiteoak, 1966; Gardner and Davies, 
1966; Gardner et al, 1967; van de Hulst, 1967; Berge 
and Seielstad, 1967; Mathewson and Nicholls, 1968]; 
and (5) the Faraday rotation of pulsar emissions [Smith, 
1968a, b; Radhakrishnan et al, 1969; Ekers et al, 1969; 


Staelin and Reifenstein , 1969; Goldstein and Meisel , 
1969; Hewish , 1970], and of other polarized galactic 
sources [Morris and Berge, 1964; Milne, 1968] , Of these 
various types of observation, only the 21-cm Zeeman 
effect and the pulsar Faraday rotation measurements 
provide a value for the magnetic field strength, and these 
refer only to some weighted mean of the line-of-sight 
component. Magnetic fields of the order of 10 pG have 
been found in some absorbing clouds; however, it seems 
likely that the field strength in such clouds is not typical 
of the interstellar magnetic field as a whole, and is in 
fact amplified as a result of the contraction of the clouds 
[Verschuur, 1969c] . The pulsar observations provide the 
most satisfactory results for our purposes since they 
refer to a region relatively close to the sun. If it is pos- 
sible to obtain reasonably good values for the rotation 
measure (RM= 0.85 }n e B • d/, with B in pG and l in 
parsecs) and the dispersion measure (DM = fn e dl) from 
the ratio RM/DM , one can obtain a value of the mean 
line-of-sight component of B, weighted with respect to 
the electron density. The measured values of B y range 

from 0.7 (iG to 3 jjlG, and it seems reasonable therefore 
to assume that the largest value is representative of the 
interstellar magnetic field strength. This value is consis- 
tent with that found using the mean value of n e 
obtained from pulsar measurements together with the 
rotation measures of extragalactic sources [Davies, 
1969]. Furthermore, a magnetic field strength of this 
order is required to account for the nonthermal galactic 
radio emission [Felten, 1966 ; Webber, 1968; Golds tein 
etal, 1970a] , although there is some uncertainty con- 
cerning the unmodulated spectrum of cosmic ray elec- 
trons. These and other arguments concerning the 
interstellar magnetic field have been reviewed by 
Burbidge [1969] . 

Interstellar atomic hydrogen can be detected by 
observations of (1) the 21-cm line in emission and 
absorption [Kerr and Westerhout, 1965; Kerr, 1968, 

1969] , including absorption of pulsar emissions [deJager 
etal, 1968; Guelin etal , 1969; Gordon etal , 1969; 
Gordon and Gordon , 1910; Hewish, 1970] ;(2) Lyman a 
absorption in the ultraviolet spectra of nearby early type 
stars [Morton, 1967; Jenkins and Morton, 1967; 
Carruthers , 1968, 1969, 1970a; Morton etal, 1969; 
Jenkins etal., 1969; Smith, 1969; Wilson and 
Boksenberg, 1969] ; (3) the absorption of soft Xrays 
from discrete sources [Gorenstein etal , 1967; Fritz 
etal, 1968; Rappaport etal, 1969; Grader etal, 

1970] ; and (4) the backscattering of solar Lyman a 
[Chambers etal, 1970; Barth , 1970a; Thomas and 
Krassa, 1971; Bertaux and Blamont, 1971]. Observa- 
tions of 21-cm radio emission suggest that the mean 
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density of neutral hydrogen in the vicinity of the sun is 
0.7 cnf 3 . The 21 -cm absorption measurements 
yield an estimate for <n jj>/ (where / is the distance to the 
source), provided one assumes a value for the spin 

temperature T s ; it is generally assumed that T s « 100° K. 
However, discrete sources such as pulsars do not show 
the absorption one would expect on the basis of 
measurements made using extended sources. This 
suggests that the interstellar medium may have a “raisin 
pudding” structure, with most of the 21 -cm absorption 
occurring in small, dense clouds with low spin tempera- 
ture [Clark, 1965; Kerr, 1969; Gordon and Gordon, 
1970] . The ultraviolet absorption measurements show 
much the same effect; that is, the column density 
obtained from Lyman a absorption lines in the spectra 
of early-type stars is usually, but not always, about 
one-tenth the value obtained from 21 -cm emission 
measurements made in the same region of the sky 
[Carruthers, 1970ft for a review of this topic] . Where 
Lyman oc and 21 -cm absorption measurements have been 
made for essentially the same path (between the sun and 
the Orion nebula) the results suggest that T s ~ 20° 
[Carruthers, 1969] . 

There seems to be a discrepancy between the 21 -cm 
emission and Lyman oc absorption measurements that 
needs to be resolved. Since the 21 -cm measurement 
refers to the entire line of sight, it is possible that the 
answer is simply that there is a low density 
(fl|j « 0.1 cnf 3 ) region near the sun, which affects the 
absorption measurements in at least some directions. In 
contrast, Jenkins etal [1969] have reported values of 
<n^ j> from Lyman a absorption measurements that 
exceed those obtained from 21 -cm emission measure- 
ments. It is possible that the 21 -cm emission process is 
not completely understood [Fischel and Stecher, 1967 ; 
Storer and Sciama, 1968] ; however, more detailed com- 
parisons are needed before we can be sure that a real 
discrepancy exists. The soft X-ray absorption measure- 
ments indicate that in the direction of the Crab nebula, 
0.3 cnf 3 [Rappaport etal, 1969; Grader etal, 
1970] , which is consistent with the 21 -cm measure- 
ments. However, in the direction of ScoXR-1, where 
Jenkins etal [1969] report ~ 1 .4-3.0 cm” 3 , the 

soft X-ray absorption measurements are somewhat 
contradictory [Fritz etal, 1968; Grader etal, 1970]. 
As far as the region in the immediate vicinity of the sun 
is concerned, recent measurements of the intensity of 
backscattered solar Lyman a can be interpreted as 
implying that <n pj> ^0.03“0. 12, depending on the model 
used [Blum and Fahr, 1970a; Thomas, 1971] . This is a 
relatively low value, but it is not inconsistent with other 


observations (all of which refer to path lengths of the 
order of 100 parsecs or more), and indeed is quite con- 
sistent with most of the ultraviolet absorption 
measurements. 

It has been conjectured that the interstellar gas might 
contain a significant fraction of molecular hydrogen 
[Gould and Salpeter , 1963; Gould et al, 1963] , and in 
cool, dense, dusty regions where the presence of a 
number of radicals has been detected from their radio 
emissions, this is to be expected [Hollenbach etal , 
1971]. Ultraviolet absorption measurements faded to 
show the presence of molecular hydrogen initially 
[Carruthers, 1967, 1970ft] ; however, recent measure- 
ments indicate that in some regions the concentrations 
of molecular and atomic hydrogen are comparable 
[Werner and Harwit, 1968; Carruthers , 1970h] . In the 
vicinity of the sun, however, where the atomic hydrogen 
density is quite low, and there is no evidence for dust, it 
is to be expected that molecular hydrogen is not an 
important constituent of the interstellar gas. 

In H II regions, the temperature of the interstellar gas 
is usually 5-10X10 3 °K, and is determined by a balance 
between heating due to photoeleetrons, and cooling due 
to recombination, collisional excitation, and free-free 
transitions [Spitzer, 1968a, b] . In H/ regions, the situa- 
tion is more complex; taking into account the effect of 
heating by low-energy cosmic rays [Hayakawa etal, 
1961; Field, 1962; Spitzer and Tomasko, 1968; 
PikeVner, 1968] and cooling by collisional excitation it 
is found that two thermally stable states are possible, 
one with neutral hydrogen densities of the order of 
0.1 cm" 3 , T ^ 10 4 °K, and ^0.1, and the other 

with neutral hydrogen densities of the order of 
10 cm -3 , r<300° K, and n e /« H «10~ 3 [Field etal, 
1969] . There are, however, some difficulties associated 
with cosmic ray heating [Werner etal, 1970] ^specially 
in view of our lack of knowledge of the galactic cosmic 
ray spectrum at low energies [Goldstein etal, 1970h; 
Gleeson and Urch, 1971ft] . Limits placed on the low 
energy cosmic ray flux by indirect observations [Habing 
and Pottasch, 1967; Greenberg, 1969; Fowler etal, 
1970; Solomon and Werner, 1970] appear to make 
cosmic ray heating ineffective. It has been suggested that 
heating by soft X rays may be more important [Silk and 
Werner, 1969; Sunyaev, 1969; Werner et al, 1970], and 
there are observational means of distinguishing between 
the two modes of heating [Silk and Brown, 1971; 
Bergeron and Souffrin, 1971; Habing and Goldsmith, 
1971] . In any case, there is evidence from 21 cm and 
other observations in support of a two-component 
model of interstellar HZ regions with properties as 
described above [Clark, 1965; Radhakrishnan and 
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Murray , 1969; Riegel and Jennings , 1969; Hjellming 
et aL, 1969; Rohlfs ; 1971 ] , and this would also be com- 
patible with the Lyman a absorption measurements 
[Jenkins, 1970]. 

It is to be expected that the interstellar gas has a 
velocity component of 20 km sec" 1 from the direction 
of the solar apex. In addition, however, there is likely to 
be a component associated with a general streaming of 
the interstellar gas relative to the stellar population, per- 
haps associated with the “density wave” structure 
described by Lin etal [1969] [see also Humphreys, 
1971; and Burton, 1971] , and a “turbulent” component 
that is evident from observations of interstellar absorp- 
tion lines [Munch, 1968] and 21 -cm observations [Kerr, 
1968, 1969; Goldstein and MacDonald, 1969] . Since the 
“streaming” and “turbulent” components of the inter- 
stellar gas velocity can easily amount to 10 km sec* 1 or 
more, we should not be surprised if the direction of 
relative motion between the sun and the gas differs 
substantially from the direction of the solar apex. 
Indeed the OGO 5 measurements of scattered solar 
Lyman a suggest that the direction of relative motion 
lies close to the ecliptic plane [Thomas and Krassa , 
1971; Bertaux and Blamont , 1971; Thomas , 1971], 
although as pointed out by Blum and Fahr [1971] this 
may result in part from the asymmetry of the solar 
Lyman a emission. 

A first estimate of the radius of the region in which 
the solar wind is supersonic can be obtained by equating 
the solar wind ram pressure to the pressure of the inter- 
stellar gas and magnetic field (see p. 616). On the basis 
of the observations described above, the radius is 
nominally of the order of 100AU if the effects of 
neutral interstellar gas, the interplanetary magnetic field, 
and cosmic rays are neglected. This is well beyond the 
orbit of Pluto, the most distant of the known planets, 
and it is probably too far to be considered as a major 
objective of the proposed Outer Planet Grand Tour 
missions (fig. 1). However, according to the analyses 
presently available, all the effects that have been 
neglected in this estimate tend to reduce the distance by 
perhaps 10 to 15 percent, and hence together they could 
give rise to a significant and favorable change in our first 
estimate. In figure 1, therefore, we have shown the solar 
wind termination as occurring somewhere in the region 
50 to 100 AU from the sun, with ~50 AU being a not 
unrealistic possibility. It is recommended that a major 
effort be made to improve our understanding of this 
problem, and to provide better values for the various 
parameters involved. 


THE INTERPLANETARY MEDIUM BEYOND 
THE ORBIT OF EARTH 

To date, in situ observations of the interplanetary 
medium have been carried out only in a very restricted 
region lying close to the ecliptic, between the orbits of 
Venus (0.7 AU heliocentric distance) and Mars (1 .5 AU). 
Extensive reviews of these observations have been given 
by Hundhausen [1968a, 1970] , Ax ford [1968] , Ness 
[1967, 1968] , Lust [1967] , Davis [1970] , and Brandt 
[1970] . For information concerning other regions we 
have had to rely on inferences drawn from less direct 
observations, notably those involving comets, radio tech- 
niques, and cosmic rays [Ax ford, 1968]. Within the 
next decade we can reasonably hope for a considerable 
extension of the region of in situ observations as a result 
of the planned missions to the vicinity of Mercury, to 
Jupiter and the outer planets, and possibly out of the 
ecliptic following a Jupiter swingby. In this section we 
consider the extrapolation of theoretical models of the 
interplanetary medium into the unexplored regions at 
great distances from the sun, neglecting external influ- 
ences (notably those due to the neutral interstellar gas). 

From the theoretical point of view the solar wind is 
probably understood rather well in a qualitative sense, 
although the complexity of the phenomenon is such that 
we are unlikely ever to be able to produce a complete, 
quantitative model. Recent reviews of theoretical work 
have been given by Parker [1965a, 1969] and Holzer 
and Ax ford [1970a] . As a result of the lack of observa- 
tions we must rely on theoretical models, even if they 
are oversimplified, to extrapolate our knowledge of the 
solar wind to heliocentric distances beyond the orbit of 
Mars, and to high heliolatitudes. The results of calcula- 
tions based on one such model are shown in figures 2 
and 3 [Leer and Ax ford, 1971]. In this model, it is 
assumed that the solar wind velocity is radial every- 
where, and that the protons are heated by hydromag- 
netic waves emitted by the sun in a region with a 
characteristic radius of 5 R e [Barnes, 1968, 1969]. 
Otherwise, it is assumed that the interplanetary magnetic 
field does not affect the dynamics apart from its effects 
on the electron thermal conductivity and on the proton 
temperature anisotropy. Whatever assumptions are made 
concerning the processes that take place near the sun (in 
r<0.5 AU), it is to be expected that at large distances 
from the sun (r> 1 AU), the predictions of any steady 
flow model are essentially as shown in figure 2 provided 
the effects of the interstellar medium are not taken into 
account. In particular the solar wind speed is very nearly 
constant and the number density varies inversely as r 2 
beyond r ^ 1 AU. 
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Figure 2 Variation of the solar wind number density n 
and radial velocity u with radial distance from the sun in 
units of solar radii a, according to the model of Leer and 
Ax ford [1971]. The dashed lines refer to the case in 
which there is no heat input [ Hartle and Sturrock, 
1968 ] . 


The interplanetary magnetic field, which is controlled 
by the solar wind, by the rotation of the sun, and by 
conditions in the solar photosphere, should be such that 
at large distances the field lines form conical 
Archimedean spirals in the manner described by Parker 
[1958, 1963] ; the spirals are determined by the inter- 
section of the surfaces 

r= |j 9 = e ° ^ 

where u is the solar wind speed, £1 is the angular velocity 
of the sun, and (r,0,0) are spherical polar coordinates. 
The angle made by the spiral relative to the radial direc- 
tion (the “garden hose” angle) is given by 


\jj = tan* 1 



( 2 ) 


The garden hose angle is shown as a function of radial 
distance and heliolatitude in figure 3(a). The compo- 
nents of the magnetic field can be written 



B Q a 2 { 2 sin 0\ 
ru J 


(3) 


where B 0 is the field strength at the corresponding point 
at the surface of the sun (r = a). The magnetic field 
strength is also shown as a function of radial distance 




Q b c 




Figure 3 Contours of various quantities in the interplanetary medium in a polar plot in 
which the abscissa lies in the ecliptic plane and the ordinate coincides with the axis of 
rotation of the sun, according to the model of Leer and Axford [1971 ] . 
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and heliolatitude in figure 3 (a). It should be noted that 
in the ecliptic plane (0 = tt/2) the direction of the mag- 
netic field becomes essentially azimuthal at large dis- 
tances from the sun (i.e., \p > 80° for r > 5 AU), and the 
field strength varies as 1/r. In contrast, in the polar 
regions (0 = 0,7r) the magnetic field lines are radial and 
the field strength varies as 1/r 2 . As a result of this 
behavior the ratio of solar wind dynamic pressure to the 
magnetic field pressure 

Smmu 2 

B 2 (4) 

tends to a constant value (j^ « 500) at large values ofr 
near the ecliptic plane, but increases as r 2 at high helio- 
latitudes. Here and elsewhere, we have taken the mean 
mass of the solar wind ions in to be 2X1CT 24 gm to 
allow for the presence of helium and heavier elements. 

Wherever the garden hose angle is small the electron 
temperature is dominated by heat conduction; hence 



where T eo is the electron temperature at some reference 
level r-r G . As the direction of the magnetic field 
becomes more azimuthal however, the effects of adia- 
batic expansion become more important, and ultimately 
T e <xr ~~* /3 . The electron temperature is shown as a func- 
tion of heliocentric distance and heliolatitude in fig- 
ure 3(6). As shown in figure 4, the electron collision fre- 
quency is everywhere comparable with, or greater than, 
the characteristic expansion rate (v exp = -(u/n)(dn/dr)), 
and hence we do not expect that the electron 
temperature should be anisotropic anywhere under 
normal conditions. In contrast the proton collision 
frequency may be less than i> exp beyond a few solar 
radii and remain so until the magnetic field direc- 
tion becomes essentially azimuthal. Thus, we can expect 
the proton temperature to be anisotropic, since when it 
approaches a constant speed the expansion of the solar 
wind takes place anisotropically (the radial dimension of 
an element of fluid remains unchanged, while the trans- 
verse dimensions expand in proportion to r). In regions 
where the magnetic field direction is essentially radial we 
should find 7^11 > T p \ while in regions where the field 
is essentially azimuthal we should find that T p ^ ^ 1/r 2 , 
T p 1 °c 1 jr so that eventually T p J- > T p ^. However, in the 
latter case the mean temperature T p = (T p ^ + 2T p ^)/3 
may decrease sufficiently rapidly for the proton collision 
frequency to become comparable to the expansion rate 



Figure 4 The characteristic expansion rate (v exp ) y the 
electron collision frequency (v ee ), and the proton colli- 
sion frequency (v pp j in the ecliptic plane as functions of 
distance from the sun in units of solar radii , according to 
the model of Leer and Ax ford [1971]. 

once more, so that collisions suppress the anisotropy and 
T p &f~ 4/3 . These effects are all evident in figures 3(c) 
and 3(d ) , which show the proton temperature aniso- 
tropy and the mean proton temperature, respectively, as 
functions of heliocentric distance and heliolatitude. 

The model from which the above description of the 
behavior of the solar wind is drawn involves several 
assumptions that are not strictly valid. The most impor- 
tant of these is that the flow is assumed to be steady and 
(near the sun) spherically symmetric. Unsteady flow 
may affect the behavior of the solar wind at great helio- 
centric distances since the interaction of fast and slow 
streams beyond the orbit of the earth can result in a 
mean flow that is considerably hotter than suggested in 
figure 3(d), with correspondingly smaller proton collision 
frequencies and hence a larger temperature anisotropy 
than suggested in figure 3(c). Departures from spherical 
symmetry that do not involve unsteady effects must 
occur since at any given time solar active regions tend to 
be concentrated in latitude; however, apart from pro- 
ducing a variation with heliolatitude of the asymptotic 
solar wind speed and perhaps slightly changing the 
temperature contours in figures 3(6), (c), and (d), it is 
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unlikely that qualitatively significant effects can arise in 
this way. 

The effects on the solar wind of the interplanetary 
magnetic field and of the rotation of the sun have been 
neglected in the model calculation apart from the con- 
trol exerted by the field on the electron thermal conduc- 
tivity and on the proton temperature anisotropy. Close 
to the sun, however, where the solar wind speed is com- 
parable with or less than the Alfven speed, the magnetic 
field is capable of inducing significant nonradial motions 
of the plasma [Weber and Davis , 1967] . Indeed, it seems 
likely that the solar wind has almost the full angular 
velocity of the sun out to perhaps 20 R Q . This rotational 
component of the solar wind velocity should decay 
inversely with heliocentric distance and become negli- 
gible at the orbit of earth and beyond. A meridional 
component of the solar wind velocity should be induced 
near the sun, causing the magnetic field to bend toward 
the solar equatorial plane as suggested by observations of 
radio sources near the sun [Hewish and Wyndham , 1963 ; 
Slee, 1966]. Accordingly, the solar wind in the polar 
regions must expand more rapidly than spherical geom- 
etry would suggest, and hence it can reach a higher 
velocity at a given heliocentric distance than it would 
achieve at the same distance in the solar equatorial 
plane. It seems possible that the observations of 
Dennison and Hewish [1967] , which suggest that the 
solar wind speed increases at high heliolatitudes, could 
be explained in this way. It is unlikely, however, that 
any of these effects lead to substantial departures from 
the description given in figures 2 and 3 of the behavior 
of the solar wind at large distances from the sun. 

It has been assumed in the previous discussion that 
the solar wind can be considered to be a perfect electri- 
cal conductor, and hence the magnetic field behaves as if 
it were frozen into the plasma. Field lines cannot 
reconnect in these circumstances, and hence adjacent 
regions of inward- and outward-directed fields are sepa- 
rated by neutral sheets that spiral away from the sun 
indefinitely in the same manner as the field lines them- 
selves. These neutral sheets can be identified with 
“sector boundaries” [ Wilcox and Ness , 1965; Wilcox, 
1968] , and their observed persistence can be taken to 
indicate that the configuration is stable against reconnec- 
tion in the sense that this would not lead to a state of 
lower energy [Axford, 1967] . As shown in figure 3(e), 
the parameter 0, defined as 

Smk(T D + T e ) 

t-sf — w 


is an increasing function of heliocentric distance in 
regions where the magnetic field lines are almost radial 
(i.e., j3 oc r l 2n ) and a decreasing function of heliocentric 
distance in regions where the field lines are almost 
azimuthal and the plasma is cooling adiabatically 
(j 3 <x r“ 4/3 ). If the stability of neutral sheets against 
reconnection is determined by 0, it seems possible that 
the sector boundaries may begin to break up somewhere 
beyond the orbit of earth where j3 is small [cfi Davis , 
1970] . This process may be aided by the divergence of 
the flow, which causes elements of plasma in the neutral 
sheet to drift apart where the magnetic field is azi- 
muthal; in contrast, near the sun where the field is 
nearly radial, the separation between elements of plasma 
in the neutral sheet stays almost constant. It is empha- 
sized that if such an instability should occur it is likely 
to be of a macroscopic rather than a microscopic nature, 
since it has been shown that quasilinear diffusion can 
easily stabilize a neutral sheet against microscopic 
tearing instabilities [Biskamp et ui. , 1970] . 

We have neglected the possible effects of plasma 
instabilities associated with anisotropies of the velocity 
distributions of solar wind particles. Perhaps the most 
important of these is the instability resulting from the 
asymmetry of the electron distribution function associ- 
ated with heat conduction [Forslund, 1970] ; this can 
reduce the effective electron thermal conductivity and 
also alter its dependence on temperature, so that the 
distribution of electron temperature could be signifi- 
cantly different from that shown in figure 3(6). Other 
instabilities may be associated with anisotropy of the 
proton temperature [ Parker , 1963; Kennel and Scarf, 
1968; Hollweg and Volk , 1970] ; however, the anisot- 
ropy is not observed to be large [Hundhausen etal , 
1967; Eviatar and Schulz, 1970] nor is it expected to 
be so (see, e.g., fig. 3(c)). Thus any instabilities that may 
result from this source might be expected to .produce 
only a rather low level of plasma turbulence. On the 
other hand, there is at least indirect evidence for impor- 
tant instabilities involving the protons and other ions, in 
that the heavier ions usually appear to have the same 
bulk speed, but a higher temperature than the protons 
[Snyder and Neugebauer, 1964; Robbins et at, 1970; 
Bame et al, 1970; Egidi et al , 1970] . Furthermore, it is 
possible that the solar wind ions are heated via the elec- 
trons, rather than directly as we have assumed, in which 
case some form of plasma turbulence must be generated 
to provide sufficiently high electron-proton energy 
exchange rates and to reduce the proton temperature 
anisotropy to the observed low values [Nishida, 1969; 
Toichi, 1971] . As far as the solar wind at great distances 
is concerned, plasma turbulence may play an important 
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role in providing an effective interaction with ions of 
interstellar origin [Holzer and Ax ford, 1971]. 

As a first approximation, the composition of the solar 
wind might be expected to be that of material having 
something like the cosmic abundance of elements in 
ionization equilibrium at the temperatures of 
1-2X1 0 6 °K that exist in the solar corona [cf Holzer 
and Axford ', 19706], The observed variability of the 
helium abundance in the solar wind [Robbins etal, 
1970; Hundhausen , 1970] and some theoretical work 
[Jokipii, 1966a; Geiss etal. , 1970a] suggest that the 
approximation might not be wholly adequate, although 
there is no indication that it is seriously wrong [Bame 
etal , 1970; Geiss et al, 19706]. As pointed out by 
Hundhausen etal [1968], the ionization state of the 
solar wind plasma is essentially “frozen” at a heliocen- 
tric distance of at most a few solar radii, where the 
electron temperature is not substantially less than it is in 
the lower corona. At great distances from the sun, radia- 
tive recombination times are very long indeed and one 
cannot expect this process to play an important role in 
determining the dynamics of the solar wind. For 
example, for an ion with charge Ze in the region beyond 
r = 1 AU where T e cc r _4/3 , the radiative recombination 
time r rec [Allen, 1963] is such that 

3X10 9 T e 3/A u\dn/dr\ 2XlQ 7 ^ 
T rec v exp ^ Z 2 n 2 ^ Z 2 

which is a large number even for the more highly ionized 
constituents (Z ^ 15). In high latitude regions where 
T e °c r -2 7 7 , r ree ^ ex *p is larger and increases with 
increasing r. Nevertheless, changes in the state of ioniza- 
tion may occur as a result of neu tralization of solar wind 
ions on interplanetary dust grains [Bandermann and 
Singer ; 1965; Banks, 1971], evaporation of the dust 
grains themselves followed by photoionization [Banks, 
private communication, 1971] , evaporation of material 
from comets and planets, and most importantly the 
penetration of interstellar gas into the inner solar sys- 
tem, as discussed in a later section (p. 623). The sugges- 
tion that the solar wind should contain a substantial 
(solar) neutral component [Akasofu, 1964; Wax etal, 
1970] encounters difficulties [Brandt and Hunten, 
1966; Cloutier, 1966’, Axford, 1969] , and would lead to 
a much larger flux of He + ions than is observed. 

INTERACTION OF THE SOLAR WIND WITH 
IONIZED INTERSTELLAR GAS AND 
MAGNETIC FIELD 

Here we consider the nature of the interaction between 
the solar wind and the local interstellar medium for the 
case where the latter is part of an HZ7 region containing 


only a negligible neutral component. Many of the essen- 
tial features of this interaction were discussed by Davis 
[1955, 1962] in his work concerning the effect of the 
solar wind cavity on cosmic rays. A more detailed anal- 
ysis was carried out by Parker [1961] , with the assump- 
tion that the effects of the interplanetary (solar) 
magnetic field can be neglected. This discussion is based 
largely on Parker’s work [1963] . 

It was pointed out by Gauser [1960] and Weymann 
[1960] that the solar wind should undergo a shock 
transition so that it can come into equilibrium with the 
interstellar medium [McCrea, 1956 ; Holzer and Axford, 
1970a] . The situation is closely analogous to the dis- 
charge from a Laval nozzle into the atmosphere 
[Liepmann and Roshko, 1957] , and to the flow of 
water over a dinner plate (see appendix). In the absence 
of a neutral component of the interstellar gas, a shock 
must always exist if the pressure exerted by the inter- 
stellar medium is finite. The suggestion made by Faus 
[1966] that the solar wind might terminate as a “free 
expansion” is incorrect; indeed, the solar wind itself can 
be regarded as a free expansion within the shock transi- 
tion. 

ej 

If the effects of the interplanetary magnetic field are 
neglected, and the solar wind is assumed to behave as a 
perfect fluid with specific heat ratio 7 = 5/3, then since 
the Mach number (M) of the supersonic solar wind is 
very large at the expected position of the shock transi- 
tion, the shock is ‘‘strong.” In this case the 
Rankine-Hugoniot relations yield 

7 - 1 7+1 22 

— P 2 ai — P l u m 


and 


K* = 


(7-1)a4 +2 

2 7 yi4-( T -l) 


7-1 

27 


( 8 ) 

( 9 ) 


Here subscripts 1 and 2 denote conditions on the 
upstream (supersonic) and downstream (subsonic) sides 
of the shock, respectively, and the subscript n indicates 
that the component of velocity normal to the shock is 
involved; p is the total pressure, and p = nm is the mass 
density of the gas. In the subsonic region beyond the 
shock, the Bernoulli equation holds along each stream- 
line: 


1 u 2 2 P = 1 P 1 
2 7-1 p 7-1 P s 


(10) 


together with 
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p7 


Ps 



( 11 ) 


where the subscript s denotes the stagnation value 
(u = 0). On combining equations (10) and (1 1) and using 
the definition M 2 - u 2 jc 2 = pu 2 /^p, the following 
equations can be obtained for p, p, c 2 , and T : 



It’ is evident from (9) that M 2 is small on the subsonic 
side of the shock transition, and if it remains small it is 
permissible to neglect terms 0(M 4 ), so that (10) reduces 
to 

1 o 

-p s u 2 +p = p s (15) 


Also from (13) we see that p = p s +0(M 2 ), so that if 
M 2 « 1 the flow can be considered as being 
incompressible. 

In general the shock is not spherical, and hence the 
flow in the subsonic region contains vorticity, and the 
stagnation pressure is not the same along all streamlines. 
However, for streamlines that lead to stagnation points 
at the interface between the solar plasma and the inter- 
stellar medium, the conditions for equilibrium require 
that p s =p 0 , where p Q is the local pressure of the inter- 
stellar gas and magnetic field. In the simple configura- 
tions of interest, these streamlines are normal to the 
shock and hence u ni =u x (the solar wind speed on the 
supersonic side of the shock). For these streamlines it 
can be deduced from equations (8), (9), and (1.0) that 



where, for 7 = 5/3, K ^ 1.13. Since n x = n e r e 2 /r s 2 i 
where r s is the radial distance to the shock and the 
subscript e denotes quantities evaluated at the orbit of 
earth (r e = 1 AU), this equation can be used to deter- 
mine r s . Note that this should be expected to be the 
minimum distance to the shock, and that r s may be 
substantially different on other streamlines. 


Let us consider the case in which the interstellar gas is 
at rest with respect to the sun and there is no interstellar 
magnetic field. The configuration is spherically symmet- 
ric, and hence (16) must apply to all streamlines if we 
interpret p Q as the pressure of the interstellar medium 
(and therefore of the solar wind) at infinity. Since the 
Mach number in the subsonic region can only decrease 
with heliocentric distance (M 2 <M 2 ^ 0.2), the approx- 
imation of incompressibility must be quite accurate. 
Hence, from mass conservation we obtain 

2 

r>r s (17) 

and thus at large distances the sun would appear to be a 
point source of incompressible fluid. 

Next let us consider the case in which there is an 
interstellar magnetic field which is uniform and has 
strength B Q at great distances from the sun. We gener- 
alize Parker's [1963] treatment of this problem by the 
inclusion of static interstellar gas, which exerts a uni- 
form pressure pg , and also by allowing for the effects of 
compressibility in the region of subsonic solar wind 
flow. It is assumed that the magnetic scalar potential in 
the external region (fig. 5) is given approximately by 



“ ~[(l/2)/ 3 /r 2 + r] B 0 cos 6 (18) 


where (r,6,0) are spherical polar coordinates with 0*0 
the axis of symmetry, and l is the radius of the diamag- 
netic sphere that would produce the same external field. 
The magnetic field lines are given by 



where r = b is the radial distance to the point where the 
field line intersects the plane 0 = 7r/2. If the interface 
between the interstellar medium and the solar wind is 
determined by (19) with b = b Q} then at large distances 
from the sun (z -r cos 0 large) the distance of the inter- 
face from the axis of symmetry is given by 

~| 1 (2 

( 20 ) 
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Figure 5 Calculated shape of the region of subsonic 
solar wind flow (shaded), for the case when the solar 
wind plasma is entirely confined by a uniform inter- 
stellar magnetic field. The solar wind is supersonic in 
r<r s , and exits in the ±z directions to infinity. The 
configuration of the interstellar magnetic field is indi- 
cated by some sample field lines in the ou ter region. The 
parameters used in this example are as follows: 
B g - 6pG , n s = 0.05 cm “ 3 , T g = 10* °K, with 
B e - 3X 10 5 G,p e - 10“ 23 gm cm~ 3 , u e = 400 km sec~ l * 

The condition of pressure balance at the stagnation line 
r = b, 0 = 7r/2, yields 



Similarly, along the interface of the exit channels, where 
M 2 -r M J; as z °°, we require for pressure equilibrium 
that 


Ps 


("'-T 


M, 


2 \-7/(7-i) 


V 


+ P e (22) 


W 


Finally, from mass conservation, and assuming that the 
flow in the exit channels becomes uniform as z 

/ y -i \ l/ 2 (y + 1 )l(y~ l ) 

2Trc 0 2 (7P s p s ) 1/2 M 00 [l+—M 00 2 J =Q 

(23) 


where Q is the total solar wind mass flux. Given B 0 , p g 
Moo, and the parameters characterizing the supersonic 
solar wind, we can solve equations (20) through (23) for 
c Q , b Q , r s , and l. As in the analogous problem of the flow 
due to a point source in a rotating fluid [Barua, 1955] , 
it is necessary to specify the exit conditions to obtain a 
unique solution. In the case of the solar wind, the appro- 
priate condition would seem to be that the fluid should 
escape freely to infinity, and hence we have taken 
Hj = 1 together with reasonable values of B 0 and p g to 
obtain the solution shown in figure 5. 

The above solutions are likely to be inadequate as 
representations of the actual interaction between the 
solar wind and the interstellar medium, since any relative 
motion between the sun and the interstellar gas must 
cause the solar wind plasma in the subsonic region to be 
blown away. This effect is illustrated in the third prob- 
lem treated by Parker [1961] : the case in which there is 
no interstellar magnetic field but a steady interstellar 
wind having a small Mach number relative to the sun. 
The solution to this problem is obtained easily if it is 
assumed that the solar wind can be replaced by a point 
source of incompressible fluid, and that the interstellar 
gas can also be considered as being incompressible. On 
combining the scalar velocity potentials for a source and 
a uniform flow, the streamlines of the resulting flow are 
found to be given by 


(z* +**)"* ' V +C ( } 

in terms of cylindrical polar coordinates (s,0,z), where C 
is a constant on each streamline, and r Q = 
^ s (P 2 ^ 2 2 /p^f /2 ) i/4 - The streamline defining the inter- 
face between the solar wind and the interstellar gas is 
given by C = - 1 if the interstellar wind flows in the posi- 
tive z direction. The stagnation point occurs at s = 0, 
z = ~r Q , and for z »r Q , the solar wind is confined to a 
circular cylinder of radius 2 r , as shown in figure 6(a). 

The major defects of this model, as far as the solar 
wind is concerned, are (1) the assumption that the inter- 
stellar gas behaves as an incompressible fluid, and (2) the 
implicit assumption that r s 2 « r Q 2 ; that is, 
u 2 /V»(p g lp 2 ) 1/2 . Although it is not easy to satisfy 
these requirements with acceptable values for the various 
quantities involved, we expect that the general nature of 
the flow remains more or less the same, as indicated in 
figure 6(a), even for the extreme case of a hypersonic 
interstellar wind. Thus, the value of r s can still be calcu- 
lated from (16), provided p Q is taken to be the stagna- 
tion pressure of the interstellar wind (p Q =p g + p g V 2 ). 
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Figure 6(a) Streamlines in a model of the interaction 
between the solar wind and the interstellar wind in 
which both fluids are treated as being incompressible 
[ Parker , 1961, 1963 ] . 

This will represent only the minimum value of r s , how- 
ever, since the cavity within which the solar wind is 
supersonic must become progressively more asymmetric 
as the Mach number of the interstellar wind increases* 
Provided we restrict our attention to the flow in the 
vicinity of the axis of symmetry in the upstream direc- 
tion, it should be possible to obtain estimates for the 
radius of curvature of the solar wind shock and for the 
positions of the interface and of a possible shock in the 
interstellar medium, using the approximate theory for a 
hypersonic shock standoff region due to Lighthill 
[Hayes and Probstein , 1966; pp. 232-254] , If the radius 
of curvature of the solar wind shock exceeds r s on the 
axis of symmetry as indicated in figure 6(6), the shock is 
not normal to the flow everywhere, and the vorticity 
induced in the subsonic region can play an important 
role in turning the flow into the wake. 

An alternative procedure using Busemann’s method 
[Hayes and Probstein , 1966; pp, 78, 156] has been 
adopted by Baranov et al. [1970] , who treat the sub- 
sonic region as a thin layer separating two hypersonic 
streams. In this case the details of the region of subsonic 
flow are ignored and one simply writes equations repre- 
senting conservation of mass and momentum for the 
layer as a whole, which can be represented as a curve of 
the form r = r s %(6). Thus, the mass flow is 

$ = Tn^PgV 2 + 27ir 2 (l - cos 0)p 2 u 2 (25) 

at a point (r,6) in the layer. The momentum equations 
are 


Pg V n P 2 U ln + 


<frv 

277 rR sin 0 


(26) 



Figure 6(b) Calculated shape of the interaction region 
of the solar and interstellar winds assuming hypersonic 
flow [Baranov et al, 1970 ] . The shape of the stream- 
lines, and the positions of shock fronts are indicated by 
the ( sketched ) dashed lines. 


— (4>v) = 27 tt sin 0(pgV n V t + p 2 tt 2 n w 2 f) (27) 
ds 6 

where subscripts n,t indicate normal and tangential com- 
ponents, respectively, v is the mean speed of the com- 
pressed gas in the layer, s is the distance measured along 
the layer from the stagnation point at 0 = 0, and 

(r 2 +r' 2 ) 3/2 

R = — • — : 

r 2 + 2r' 2 -rr" 

is its radius of curvature. Note that the equation for 
momentum conservation in the direction of the normal 
contains a term representing the centrifugal force acting 
on the plasma in the layer, but this does not affect the 
distance to the stagnation point (r = r s ) since 4>w== 0 at 
0 -0. With some manipulation v and <$> can be elimi- 
nated from these equations, leaving a third-order, non- 
linear differential equation for £(0) that must be solved 
numerically (see fig. 6(6)). The thickness of the layer, 
which can be obtained using the Lighthill method, 
would be expected to be of the order of 30 percent of 
the radius of curvature of the layer at the stagnation 
point (R = 5r s /3). 

The closely analogous problem of the interaction 
between the solar wind and the ionized head of a comet 
has been treated as being equivalent to a source of 
incompressible fluid in a supersonic flow by Ioffe 
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[19 66a f b ; 1968] . Other aspects of the case in which the 
interstellar wind is supersonic have been considered by 
Dokuchaev [1964], Brengauz [1969], and Sakashita 
[1970], who note that the shock waves produced by 
randomly moving stars surrounded by extensive stellar 
wind regions could be a significant source of energy for 
the interstellar gas. It should be noted, however, that 
provided the stellar winds are supersonic near the stars, 
there can be no drag whatsoever on the stars themselves. 
The momentum imparted to the interstellar gas is taken 
up by the outer subsonic part of the stellar wind region 
as it is turned to form a comet-like wake behind the star. 
In fact, this must be essentially correct even if the stellar 
wind is subsonic everywhere, provided the dimensions of 
the star are small compared with those of the associated 
stellar wind region. 

Although we have not been able to construct a com- 
pletely satisfactory model of the interaction between the 
solar wind and an ionized interstellar medium, three 
significant results have been established: (1) in the 
presence of an interstellar wind the region of subsonic 
solar wind flow should form a comet-like tail in the 
wake of the sun; (2) the cavity within which the solar 
wind is supersonic can be expected to be asymmetric; 
(3) the minimum radius of the cavity is determined by 
the maximum pressure exerted by the interstellar gas 
and magnetic field, since the point of maximum external 
pressure must be a stagnation point of the (subsonic) 
interior flow. Thus, we can calculate the minimum 
radius of the solar wind shock transition from equa- 
tion (16) by taking p Q to be the maximum external 
pressure: 


n e mu x 






rtgilkTg + mV 2 ) 


(28) 


where a is a factor that allows for the enhancement of 
the interstellar magnetic field due to the presence of the 
solar wind (cf eq. (21)). The minimum r s according to 
(28), for a range of possible values of B 0 and rig , and 
fixed values of the remaining parameters, can be read 
directly from figure 7. It should be noted that in the 
absence of any ionized interstellar gas, r s & 140 AU if 
we assume that B 0 ^ 3X1 CF 6 gauss as suggested by the 
observations described earlier. If we assume that 
«g.^0.05 cm" 3 then with 20 km sec" 1 our estimate 
forr s is reduced to 100 AU. For the transition to occur 
in the vicinity of the orbit of Jupiter as suggested by 
Lanzerotti and Schulz [1969] , or even closer to the sun 
as suggested by Brandt and Michie [1962] , quite 



Figure 7 Minimum distance to the solar wind shock 
transition (r s ), as a function of interstellar magnetic field 
strength for various values of the interstellar electron 
density. The curves have been calculated from equa- 
tion (28) assuming a- 2.25, Tg=10*°K, and 
V= 20 km se<T l . 

implausible values of the parameters involved in equa- 
tion (28) are required, or else an entirely different model 
should be considered, as discussed later (p. 623). 

In the previous discussion we have ignored the effects 
of the interplanetary magnetic field on the flow of the 
solar wind. This is a reasonable procedure in the super- 
sonic region in the absence of effects associated with 
neutral interstellar gas, since > 500, and in any case 
the Maxwell stresses are very nearly in self equilibrium. 
However, in the subsonic region it is easily seen that 
inconsistencies arise if the effects of the magnetic field 
are neglected and the flow is treated as if it were incom- 
pressible. For example, if we assume that the flow is 
spherically symmetric and radial, and that the magnetic 
field is azimuthal, then since 

ruB = constant (29) 

we see that with u cc r~ 2 (eq. (17)), B& r and hence the 
magnetic field strength increases indefinitely. Although 
the flow is likely to be nonradial and the dimensions of 
the subsonic region are restricted as indicated in fig- 
ures 6(a) and (j b ), the tendency for the magnetic field 
strength to increase is unavoidable. Indeed, a very similar 
situation is found in the earth’s magnetosheath where 
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the interplanetary magnetic field strength near the solar 
wind stagnation point is often found to be comparable 
to that of the geomagnetic field just within the magneto- 
pause [Cahill and Amazeen, 1963] . 

The manner in which the magnetic field can be 
enhanced in the subsonic region, so that it ultimately 
dominates the plasma behavior, is demonstrated rather 
well by the case of steady, radial flow with spherical 
divergence, in the presence of a purely azimuthal mag- 
netic field [Cranfill, 1971]. Let us assume that the 
region of subsonic flow begins at. a shock transition 
occurring at r = r s : the Rankine-Hugoniot equations can 
be written in normalized form 


P2 U 2 = 1 (30) 



B 2 U 2 = l (33) 

where p =p/p t , m=u/w 1( B = B/Bi , p = p/p l u t 2 , and 
subscripts 1 ,2 denote upstream and downstream condi- 
tions, respectively. The equations of motion in the sub- 
sonic region can be integrated to yield 


~pur 2 = 1 


(34) 


2 

(35) 

1 )p & v 2 

% 

p _'P 2 


(36) 

V " Jj 


riiB = l 


(37) 


where r-r/r s . For a given value of j3 v , equations (30) 
through (33) can be solved for ~P 2 and p 2 , and equa- 
tions (34) through (37) can in turn be solved for u, B, p , 
and p as functions of r. It should be noted that the 
normalization results in the elimination of r s from the 
equations, leaving j3j, as the only parameter. The solution 
for a value of /3y appropriate to the ecliptic plane is 
shown in figure 8. It is evident that 7 must be large to 



Figure 8 Solution of equations (34)-(37) with initial 
conditions determined by equations (30f(33) and with 
I % = 500. Note that the presence of the magnetic field 

has caused the total pressure p T to decrease to zero as 

— 1 (o) 

r-»oo rather than remaining constant (pj )as it does in 

the field-free case [from Cranfill , 197 i] . 


make a substantial change in the total pressure 
Pj 2 ~p +pu 2 +B 2 lfi v , although the changes are not 
insignificant for moderate values of T. At 7=2 for 
example, p j is reduced by about 12 percent from the 
value p j obtained in the absence of magnetic field. 

In a more realistic model with nonradial flow it 
should be expected that the effectiveness of the mag- 
netic field in reducing the total pressure beyond the 
shock transition is more marked than in the model with 
strictly radial flow. Lateral flow along the magnetic field 
lines permits the gas pressure to decrease and be com- 
pensated by a corresponding increase in the pressure of 
the magnetic field. It is easily shown that the reduction 
in total pressure A p across a layer of thickness d and 
radius of curvature R is given approximately by 

where {B 2 > is the mean square magnetic field strength in 
the layer. This result applies equally well to the earth’s 
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(cm 2 sec ster ) 


magnetosheath, and one might expect the values of 
&pj p s to be comparable in the two cases. Observations 
made in the magnetosheath suggest that Ap/p s could 
easily be 0.1 ~0.2, and hence we can reduce the minimum 
distance to the solar wind shock transition from our 
previous estimate of 100 AU to perhaps 80 to 90 AU. 

We have assumed that the interplanetary magnetic 
field can be considered as being frozen into the solar 
wind plasma in the above discussion. However, magnetic 
field directional discontinuities must exist within the 
solar wind plasma (sector boundaries), and at the inter- 
face with the region of interstellar plasma and magnetic 
field. Throughout most of the region of subsonic solar 
wind flow, it might be expected that sector boundaries 
are stable against reconnection since |3 tends to be 0(1) 
or greater. However, there is no obvious reason why the 
interface with the interstellar region should be stable 
against reconnection, and hence we can expect that the 
interstellar and interplanetary (solar) magnetic fields are 
connected in much the same manner as the inter- 
planetary magnetic field connects with the geomagnetic 
field [Dungey, 1961 ; Levy etal, 1964]. 

It would be of great interest if the shell of hot plasma 



formed by the region of subsonic solar wind flow could 
be observed from a distance; however, this does not 
appear to be feasible. The suggestion that free-free 
transitions in this region might produce a detectable flux 
of soft X rays with a roughly isotropic distribution was 
originally made by Reiffel [1960] . The spectrum of the 
observed isotropic X-ray background is shown in fig- 
ure 9, together with calculated spectra based on the 
(favorable) assumptions that the electron and proton 
temperatures are equal in the shell of hot solar wind 
plasma, and that the thickness of the shell is 4 r s . The 
contribution from free-bound transitions associated with 
the ultimate recombination of all stellar wind plasma 
originating from stars within 330 parsecs of the sun is 
also shown. It is evident that neither the solar wind nor 
stellar winds in general can be considered to be likely 
sources for the observed X-ray flux [Bowyer and Field , 
1969; CranfiB, 1971]. 

The shock transition that terminates the supersonic 
solar wind may be collision free, in which case its 
thickness is no more than a few proton gyroradii (about 
10 10 cm for 6 = 7r/2, r s — 100 AU). However, it is also 
possible that the transition is as wide as the system will 


Figure 9 A comparison of observations of the diffuse 
soft X-ray flux (dashed line), with estimates of the X-ray 
flux resulting from (a) free-free emission from the region 
of subsonic solar wind flow (straight solid lines), and 
(b) recombinations of stellar wind ions (bar graph ). The 
thickness of the region of subsonic solar wind flow is 
taken to be 4r s , with T e = 184 eV t and n- 4n e r e 2 /r s 2 
with n e = 5 crn 3 . Even though we have assumed a rather 
large thickness for the subsonic region , it is evident that 
for reasonable values of r s (i.e., r s > 50 AU) the free-free 
emission fails to account for the observed flux by a 
factor of the order of 10 3 . The bar graph represents 
0.05 keV averages of the free-bound emission resulting 
from the recombination of heavy ions that have been 
injected into the interstellar medium by stellar winds . It 
is assumed that a steady state has been reached with the 
recombination rate equal to the injection rate , and that 
all stars within 330 parsecs of the sun contribute. The 
star density is taken to be 0.08 parsec 3 , and all stellar 
winds are assumed to be similar to the solar wind in 
composition and total flux. Ground state recombina- 
tions of all appropriate ion states of He, C, N, O, Ne, Mg, 
Si, and Fe have been included. It is evident that the 
recombination radiation is also insufficient to account 
for the observed soft X-ray flux [Cr an fill, 1971). 
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permit [, Schindler , private communication] , in which 
case the shock structure would be determined by 
proton-proton collisions. The thickness of such a shock 
is presumably of the order of X cos 2 \jj ~ X(m/t£ 2 sin 6) 2 
where X is the mean free path. Since at large distances 
X/u < l/v ex p » r/2u (fig. 4), it is evident that the thick- 
ness must be less than ll(2r s sin 2 0) AU. Thus, in the 
ecliptic plane (0 = tt/2), we expect that the shock thick- 
ness is less than 10 1 1 cm, even if the shock is collision 
dominated. 


NEUTRAL INTERSTELLAR GAS IN THE 
VICINITY OF THE SUN 

If there were no relative motion between the sun and the 
surrounding interstellar gas, and if the complications 
introduced by the solar wind could be ignored, the sun 
could produce a quite substantial H II region if the 
density of the ionized gas is not too large [Newkirk etal 
1960; Brandt , 196%; Lenchek , 1964; Williams , 1965] . 
However, if there is relative motion, the situation is 
drastically altered, and even with quite modest speeds of 
the order of 20 km sec” 1 it can be shown that the 
neutral interstellar gas can penetrate almost unattenu- 
ated into the inner solar system within the orbit of 
Jupiter [Blum and Fahr, 1969, 1970a; Holzer, 1970; 
Holzer and Axford , 1971; Semar, 1970] . In this section, 
we assume that the interstellar medium in the vicinity of 
the solar system is part of an H/ region, and examine the 
manner in which the neutral gas is affected by the solar 
wind and by ionizing radiation from the sun. 

Let us consider first the very simple model in which 
the nearby interstellar gas is stationary with respect to 
the sun, with the solar H/7 region assumed to be fully 
ionized and having the same pressure as the surrounding 
El region. On equating the flux of ionizing photons 
emitted by the sun Q to the total recombination rate 
we find that the radius of the H/7 region rjj is given by 
Stromgren’s result: 


r II = 



1/3 


(39) 


where a « 3X10* - 107 F is the coefficient for recom- 
bination to all levels [Alien, 1963, p. 90] . The condition 
of pressure equilibrium yields 


2 n n T u ~ n i T i < 4 °) 


2 = 4X1 0 3 7 sec” 1 ; the points of intersection of the two 
sets of lines give the values of rjj for the various com- 
binations of values of Tj and t ii- The thickness of the 
ionization front separating the H7 and Eli regions is of 
the order of A = \ jojjij, where « 6.3X10” 18 cm 2 is 
the photoionization cross section for hydrogen atoms at 
the Lyman limit. In cases where rjj< A & 5X 10 4 AU 
our assumption that the Eli region is well defined is 
invalid and a more detailed analysis must be carried out 
[Gould etal. , 1963; Lenchek, 1964; Williams, 1965]. 
The thickness of the ionization front is somewhat 
greater than A since the absorption coefficient decreases 
approximately as v 3 beyond the Lyman limit, and con- 
sequently the radiation reaching the outer part of the 
Eli region is such that the effective absorption coeffici- 
ent is much less than o Nevertheless, the value of 
calculated from equation (39) has some significance as a 
“characteristic 55 dimension, even in circumstances where 
the Eli region is not sharply bounded. 

In principle, the solar H/7 region could contribute to 
the low frequency cutoff in the spectrum of galactic 
nonthermal radio noise, as suggested by Lenchek 
[1964] . However, the density required to produce sig- 
nificant absorption at frequencies of the order of 1 MHz 
is quite large cm" 3 ) even if the electron 

temperature is as low as 10 3 ° K. 

In the absence of any interstellar magnetic field, there 
would be little point in considering a static H77 region of 
the type described above, since the solar wind would 
push the interstellar gas to such a great distance that 
ultimately the H/7 region would either overlap those of 
other stars or be limited by recombination of the solar 
wind plasma. In the latter case, the solar wind would 
terminate at a recombination front [Newman and 
Axford, 1968] . The interstellar magnetic field should 
not be expected to have any important direct effect on a 
static, approximately spherical solar Eli region; how- 
ever, it would effectively shield the Eli region from the 
solar wind, which would be channeled harmlessly away 
in the manner indicated in figure 4. 

The chief defect of these models of a solar Eli region 
is the requirement that there be no relative velocity 
between the sun and the local interstellar gas, when in 
fact a relative speed of the order of 20 km sec -1 might 
be considered more appropriate. An atom of interstellar 
gas, approaching the sun with speed V - 20 km sec" 1 , 
has a mean free path against loss by ionization given by 


Figure 10 gives these relationships for various values of 
Tj and Tjj for the case nj = 0.2 cm" 3 , with 



(41) 
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Figure 10 The values of rjj and nj (solid lines) are plotted for various values of T^and 
7 'i/Th according to equations (39) and (40). The intersection of these lines defines repre- 
sentative allowed sets of parameters describing model solar YUI regions. The ionization 
front thickness A is shown for two values of Tj/Tjj (dashed lines). As an example , for 
T n = 3X10 3 °K, and Tj= 30° K, we see that r H ** 7X10 4 AU f and A » i<9 4 AU. 


where the 1 3j are rate coefficients for all possible ioniza- 
tion processes (table 1). In general, ~ P 0 (r e 2 /r 2 ), where 
is a constant for each atomic species. It is to be 
expected therefore that a substantial fraction of inter- 
stellar atoms approaching the sun can penetrate to 
within a heliocentric distance of the order of 
r i = (fi Q r e 2 /V). From table 1 we see that interstellar 
hydrogen atoms can easily penetrate to within the orbit 
of Jupiter, while interstellar helium atoms can penetrate 
to the orbit of earth. Clearly then, we must expect to be 
involved in quite a different situation than that of a 
static, spherical UII region that effectively excludes 
neutral interstellar gas from the region occupied by the 
solar wind. 


The flow of neutral interstellar gas past the sun has 
been studied by Fahr [1968a] , Blum and Fahr [1969 , 
1970&] , Holzer [1970], Holzer and Axford [1971], 
Tinsley [1971] , and Semar [1970] . The following dis- 
cussion is based on Holzer’ s treatment of the problem 
(essentially the same as that of Blum and Fahr), with a 
slight generalization to allow for solar radiation pressure. 
As a first approximation, solar radiation pressure can be 
assumed to produce a radial outward force that varies 
inversely with the square of the distance from the sun. 
Thus it can be taken into account by introducing an 
“effective” gravitational constant (l-ju)G, where ju is the 
ratio of the force on an atom resulting from radiation 
pressure to the force of gravity. For hydrogen, the large 
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Table 1 Loss coefficients at r = rj & = & Q ), and penetration distances for 
V= 20 km sec” 1 (r^j. The solar ultraviolet flux has been taken from Hinteregger [ 1970] 


Species 

Reaction 

a(cm) 2 

Reference 

Rate 

(lcrV 1 ) 

r i(AU) 


H + hi> -»■ H + + e 


Banks and Kockarts [1972] 

1.5 


H 

H + H + H + + H 

2X1CT 15 

Fite et al [1962] 

4 

4 


H + He +2 -► H + + He + 

4X1CT 16 

Fite et al. [1962] 

0.036 



He + hr> He + + e 


Lowry et al [1965]; Cairns and 

0.62 





Samson [1965] 



He 

He + H + -*■ He + + H 

1.5X1CT 1 7 

Afrosimov et al [1969] 

0.03 

0.5 


He + He ++ -► 2He + 

2xic r 17 

Hertel and Koski [1964] 

0.0018 


C 

C + hv -*■ C + + e 


McGuire [1968] 

40 

30 

N 

N + hv -> N* + e 

1.8X10" 15 

Henry [1968] 

1.7 

A 

N + H + N + + H 

Stebbings etal [1960] 

3.6 

4 

0 

0 + hi> -+ 0 + + e 


Huffman [1969] 

2.4 

o 

0 + H + -»• 0 + + H 

8X 10" 1 6 

Fite et al [1962] 

1.6 

3 


Ne + hv -*■ Ne + + e 


Ederer and Tomboulian [1964] , 

1.9 


Ne 

Ne + H + -> Ne + + H 

10" 16 

Samson [1965] 
Afrosimov etal [1969] 

0.2 

1.6 

Si 

Si + hp -* Si + + e 


McGuire [1968] 

60 

45 


Ar + hv -*■ Ar + + e 


Cairns and Samson [1965] ; 

2.8 


Ar 



Rustgi [1 964] ; Madden et al 
[1969] 


3.6 


Ar + H + -► Ar + + H 

10" 15 

Koopman [1967] 

2 


Fe 

Fe + ht» -*■ Fe + + e 
Fe + H + -» Fe + + H 

CIO" 16 

McGuire [1968] 

Lee and Hasted [1965] 

2.4 

1.8 


solar Lyman a flux results in p being of order unity, and 
hence the effects of radiation pressure are very impor- 
tant [Wilson, 1960; Brandt, 1961; Tinsley, 1971] , which 
is contrary to the conclusion of Fahr [1968a] . However, 
for other species p « 1 , and hence solar gravity deter- 
mines their dynamic behavior. In a more accurate treat- 
ment it would be necessary to take into account the 
variation of p with position resulting from Doppler shifts 
and the reduction in the intensity of the solar radiation 
resulting from scattering. It should also be noted that 
the solar Lyman a flux varies with solar activity espe- 
cially at the center of the line where 100 percent varia- 
tions can occur [Blamont and Madfar, 1971] , 

If we neglect the thermal velocities of interstellar 
atoms (<1 km sec* 1 if Tj~ 100° K) in comparison with 
the bulk velocity of the gas relative to the sun 
(V 20 km sec* 1 ), then the trajectory of every atom lies 
in the plane determined by its velocity vector at infinity 
and the sun. Since the mean free path for atom-atom 
collisions is very large ( — 10 1 10 1 7 cm), and the most 


important collisions with solar wind particles are 
ionizing, it is appropriate to make use of a free particle 
treatment allowing for losses associated with ionization, 
with loss coefficients as given in table 1 . Taking polar 
coordinates (rfi f ) in any such plane, it can readily be 
shown that the trajectories are hyperbolae with the sun 
as focus, given by 

— = 1 + B sin(0' + a) (42) 

r 

with 

= (43) 

M0 = 7 (44) 
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where 0' - 0 is antiparallel to the incident flow velocity 
V, p is the angular momentum of the atom about the 
sun, Up uq are velocity components in the r f 0 f direc- 
tions, M is the mass of the sun, and 


r 

(45) 

1/2 


/ F 2 /! 2 ^ 

(46) 

( ,+ ) 

1 

sin a = --- 

(47) 


In general, any point in the plane (r,d’) is a point of 
intersection of two trajectories having angular momenta 
given by 

P 1 .P 2 -\v 



Figure 1 1(b) Same as (a), but with the effects of radia- 
tion pressure exceeding gravity. A parabolic forbidden 
region (hatched) is produced with the sun at the focus. 


( AGM\ 

— ) Ox-1) (49) 


(r sin 0' ± 


r 2 sin 2 0' + 


4r(I - cos#') 


(48) 

where C “ V 2 1(1 — p)GM, and the plus sign is associated 
with p x and the minus sign with p 2 . As indicated in 
figure 11, p t > 0, p 2 < 0 if p < l,andp 1 >0 ,p 2 > 0 if 
p > 1 . In the case ju < 1 , some particles may be lost by 
striking the sun. In the case p > 1, the stream of parti- 
cles cannot enter the region defined by 


since the angular momenta given by equation (48) 
become imaginary. In the case p = 1 , the trajectories are 
straight lines parallel to 0' = 0 and p 2 = 0. 

The number density n of a particular species in inter- 
planetary space can be calculated from the equation of 
continuity, with the trajectories (42) being treated as 
streamlines. We use an orthogonal coordinate system 
(£,^,0) such that e ^ = u/w, and 0 is the azimuthal coordi- 
nate relative to the axis 0' = 0. Thus 



Figure 11(a) Examples of intersecting particle trajec- 
tories for a case of a cold interstellar wind approaching 
the sun (at the origin) from the right of the diagram at 
20 km sec 1 , and with radiation pressure absent (p - 0). 
Note that all trajectories cross the axis of symmetry on 
the downwind side of the sun. 


< 5 °> 

with ft as given in table 1 , and 

f 1 d [ ln(nuh = - f 1 (~W d% (51) 

OO OO 

If we take = p/V, then it can be shown that 


h\ y = 


\2V-r sin 0'\ [r 2 sin 2 6' + 4r(l-cos 6')/C ] 1/2 


1*1(1 + 2/rC) 


, 1/2 


|*|(1 + 2 IrC) 


1/2 


Also, since the displacement d% is along a trajectory 


(52) 


, , fc _ urdd ' 


(53) 
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Hence on putting h ^ = r|sin 0' I , and taking j3 = ft/ e 2 /r 2 , 
we obtain from (51) 



exp(lV<? 2 0j/lp/l) 

K*r sin 6 [r 2 sin 2 0 + 4r(l-cos 0)/C] 1 /2 


/ 


(54) 


where 0 is the polar angle from the axis of symmetry 
(0 < 0 < 7r and Isin 0' I = sin 0), N is the number density 
of the incident stream at infinity, and Oj ~ 6 if py >0, 
0y = (2 tT“ 0) if pj < 0. If either of the trajectories that 


intersect at (rfi') strikes the sun before reaching the 
point, its contribution to n must be discounted. Note 
that if p •-* 1, then u = F, = Fr sin 0', p 2 = 0, 

and (54) reduces to the simple result 

n(r,Q) = N exp ( ) (55) 

\Frsin0 / 

Equal number density contours for hydrogen and 
helium have been calculated withp = 0, 0.8, 1.0, and 1 .2 
for the case of hydrogen (figs. 12(a)-(d)) and withp = 0 



(c) p - 1.0 (d)p= 1.2 


Figure 12 Contours of equal density of neutral interstellar hydrogen. The sun is at the 
center of each diagram and the orbit of earth is shown as a dashed line. The plane of the 
diagram is the plane containing the velocity vector of the interstellar wind and the sun. 
The three-dimensional contours are of course axisymme trie. Note the high density on the 
axis of symmetry (8 = tx) for p<l, and the paraboloidal void which appears when p 
exceeds 1. 
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for helium (fig. 13). Note that n(r ; 0) is finite on 0 = 0 in 
each case, but if /x < 1 , n(r,0) 00 as 0 it as shown in a 

slightly different context by Danby and Camm [1957] . 
This singularity disappears if one relaxes the assumption 
that the incident stream of atoms has no thermal veloc- 
ities [Danby and Bray , 1 967] . For 1, n(r,n) = 0, and 
for p > 1 an empty region develops in the wake of the 
sun. In the upstream direction the interstellar neutral gas 
penetrates to a heliocentric distance of order 
r i = $o r e IV* as argued previously, and hence there are 
likely to be observable effects that permit us to deter- 
mine N(E),N(Ee), andV. 

Fahr [1968a,h; 1969] has examined the possibility 
that the local fluxes of neutral interstellar gas and of 
neutral hydrogen produced in the solar wind by charge 
exchange should be sufficient to cause changes in the 
structure of the upper atmosphere. However, on taking 
reasonable values for the various parameters concerned, 
one finds that it is unlikely that observable effects can 
be produced. Holzer and Axford [1971] have argued 
that accretion of interstellar gas could be an important 



Figure 13 Contours of equal number density of inter- 
stellar helium in the plane containing the sun and the 
velocity vector of the interstellar gas. The orbit of the 
earth is shown as a dashed line . 

item in the terrestrial budget of 3 He but not of other 
gases unless conditions were vastly different in the past. 
It is evident from the calculations described above that 
the number density of neutral helium near the earth is 
only slightly less (^ 25 percent) than it is in the interstellar 
gas. If we assume that [ 4 He] :[ 3 He] 3000:1 
[Cameron, 1968] , it is estimated that the average flux of 
interstellar 3 He atoms accreted by the earth is approxi- 
mately 50 NQT) atoms cm“ 2 sec - 1 . According to 


Johnson and Axford [1969] , the average rate of loss of 
atmospheric 3 He by evaporation is approximately 
7 atoms cnT 2 sec -1 ; hence, allowing for contribu- 
tions from other sources, we find that Af(H) ** 
0.05—0.1 cm -3 . At certain times of the year, the 
relative velocity between the earth and the interstellar 
gas can be as high as 80 km sec - 1 . This is sufficient for 
heavier atoms to be trapped in a target foil, thus raising 
the possibility that the isotopic composition of the inter- 
stellar gas (especially neon and argon) might be mea- 
sured directly [Geiss, private communication] . 

It might also be possible to estimate the density of 
the neutral interstellar gas in the vicinity of the solar 
system from observations of the state of ionization of 
the solar wind. For example, charge exchange between 
solar wind protons and interstellar hydrogen atoms pro- 
duces a neutral hydrogen component of the solar wind; 
similarly, photoionization of interstellar helium atoms, 
together with other less important processes, results in 
the appearance of He* ions in the solar wind (fig. 14). 
The hydrogen flux in the solar wind is estimated to be 
less than 10 -4 of the proton flux if 7V(H) = 0.1 cm" 3 , 
and is probably too small to be easily measured. The 
expected flux of He* ions in the solar wind is shown in 
figure 14(d) for the case Af(H) = 1 cm -3 . There have 
been reports of a singly ionized helium component of 
the solar wind [Bame etal, 1968; Wolfe etal, 1966], 
with a flux of the order of 10“ 3 times that of He* 2 ions. 
On the basis of these results, Holzer and Axford [1971] 
have estimated that Af(H)« 0.1 cm -3 , which is consis- 
tent with the estimate based on the terrestrial 3 He bud- 
get. More recent observations by Bame et al [1970] 
showed no evidence for the presence of 4 He* ions in the 
solar wind possibly because they were masked by other 
ions with the same charge-mass ratio (Si* 7 , S* 8 , A* 9 , 
Fe* 14 ) or because they had a bulk flow velocity and 
temperature significantly different from those of other 
ions and could not be easily identified in an energy -per- 
unit-charge spectrum. There is clearly a need for further 
observations using a detector which can identify 4 He* 
ions unambiguously [Ogilvie etal , 1968; Ogilvie and 
Wilkerson, 1969] . 

The interstellar gas can be detected directly by obser- 
vations of resonant scattering of solar photons, notably 
Lyman a (XI 2 16 H/) and X584 He/ [Morton and Purcell , 
1962; Kurt , 1965, 1967; Kurt and Dostovalov , 1968; 
Kurt and Syunyaev, 1968 ; Young et al, 1 968 ; Chambers 
etal , 1970; Mange and Meier ; 1970; Thomas and 
Krassa, 1971; Bertaux and Blamont, 1971; Wallace , 
1969; Tinsley , 1969; Reay and Ring ; 1969; Barth, 
1970 a,b\ Barth etal, 1967, 1968; Metzger and Clark, 
1970; Byram et al, 1961; Bowyer et al, 1968 \Johnson 
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<b) ! 


+ (He + )/+(H*) 


3 5 AU 


(d) 



Figure 14 Contours of equal flux of hydrogen in the solar wind resulting from charge 
exchange with interstellar hydrogen , relative to the solar wind proton flux (assumed to be 
2X10 8 cm~ 2 sec~ l at 1 All). The plane of the diagram is defined by the velocity vector 
of the interstellar gas and the sun. There is an enhancement of solar wind hydrogen along 
the downwind axis of symmetry in (a) and (b) where p = 0.0 and 0.8 , respectively , but 
not in (c) where p = 1. 

The flux of He + ions in the solar wind resulting from ionization of the interstellar gas, 
relative to the solar wind proton flux , is shown in (d). The interstellar gas density has 
been assumed to be 1 cm~ 3 in these calculations , with N(Y{) = 0.92 cm~ 3 , and N(He) = 
0.08 cm” 3 . For other values of N(H) and N(He) the figures should be rescaled 
appropriately. 


etal. , 1971; Ogawa and Tohmatsu , 1971]. From such 
observations, we can in principle deduce the composi- 
tion and density of neutral interstellar gas flowing past 
the solar system, the velocity vector, and also the tem- 
perature. Accordingly, they are of great significance for 
our understanding of the behavior of the solar wind at 
great distances from the sun. Interpretations of these 
observations in terms of various models of the inter- 
action between the solar wind and the interstellar 
medium have been given by Patterson etal [1963], 
Brandt [19646, 1970] , Hundhausen [19686] , Kurt and 


Germogenova [1967] , Blum and Fahr [1969, 19706,c, 
1971] , Tohmatsu [1970] , Bowers [1970], Tinsley 
[1971], Holzer [1971], Thomas [1971], and Wallis 
[1971]. 

Let us consider a simple model for the scattering of 
solar radiation by neutral atoms in the interplanetary 
medium. The flux J(r,.fyp) of solar photons in the fre- 
quency range (v 9 v + dv) satisfies the equation 

\ ~ (r 2 J) =~anJ (56) 

r dr 
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where <x(r,d; v) is the absorption coefficient, and we have 
neglected multiply scattered photons. If r, r are the 
position vectors of a point in space relative to the sun 
and earth, respectively, then the intensity of singly 
scattered photons in a fixed direction f is given by 

f(v)=~ f a(r ; v)n(t)J (r\v)dr (57) 

4?r * 0 

where the integration is carried out along the line of 
sight from the earth, v refers to the frequency before 
scattering takes place, and absorption of the scattered 
photons is neglected. The total intensity of singly 
scattered photons from a given solar emission line is 
therefore 

/ = 47tX1(T 6 r° f{v)dv (Rayleighs) (58) 

•'-oo 

It does not matter that the frequency of the scattered 
photons need not be the same as that of the absorbed 
photons (in the observer’s frame) since observations are 
usually carried out with broadband detectors. However 
it might be important in some cases to allow for Doppler 
shifts in calculating the absorption coefficient if the 
solar emission line is not broad and featureless. In the 
case of Lyman a for example (see fig. 15), the intensity 



Figure 1 5 Solar Lyman a profiles for two regions on 
the sun obtained by Meier and Prinz [1970] . The central 
core is produced by the earth's hydrogen corona , and 
the dashed line indicates the expected unattenuated 
intensity . 


at the center of the line is only 70 percent of the maxi- 
mum intensity in quiet solar conditions, and the maxi- 
ma occur at ±0.25 A from the center (corresponding to 
Doppler shifts such that \u r I ^ 60 km sec" 1 ). If it is con- 
sidered that the above analysis is unnecessarily compli- 
cated in view of the fact that the solar emission line is 
quite variable, and that multiple scattering has been 
neglected, a simplification can be achieved by assuming 
that the medium is optically thin. In this case the right 
hand side of equation (56) is equated to zero, so that 
instead of 

J(r;v) = ■— J(r e ;v)exp ' a(r;iO«(r)cfr| (59) 

with the integration being carried out along the radius 
vector from the sun, we neglect the exponential term 
and simply use 

r 2 

J{r,v) = ~J(r e> v) (60) 

r 2 

in (57). In any case, these calculations will usually tend 
to overestimate the intensity of the scattered light for 
given N , as a result of the neglect of multiple scattering. 

Examples of calculated total intensity profiles of 
scattered solar radiation from neutral hydrogen 
■(XI 2 16 H/) are shown in figures 1 6(a) and (b) and from 
neutral helium (\584Hei) in figures 16(c) and (d). In 
figures 1 6(a) and (c) the earth is shown situated on the 
upstream axis of symmetry (a e = 0), and in figures 16(b) 
and (d) it is shown situated at 90° from the axis of 
symmetry (a e = 90°, 270°). We have not attempted to 
obtain an exact fit between the models and the observa- 
tions (fig. 17), but instead have simply taken 
N( H) = 0.1 cm" 3 , jV(He) = 0.008 cm” 3 , and 
V = 20 km sec" 1 , with the interstellar wind vector in the 
plane of the ecliptic; this choice provides a reasonable 
match to the observations if /x « 1 . The direction of the 
interstellar wind is determined by the direction of maxi- 
mum Lyman a intensity, although as pointed out by 
Blum and Fahr [1971], some distortion may be pro- 
duced as a result of asymmetries in the solar Lyman a 
emission. The velocity of the interstellar wind is deter- 
mined by the parallax effect shown in figure 16(b) (see 
also figs. 17(a)— (c)). The value of AT(H) is determined 
from the maximum intensity of scattered Lyman a, 
although it should be noted that there is some ambiguity 
since it is also necessary to account for the observed 
minimum intensity of ~250 R [ Thomas , 1971]. Note 
that for small values of p the peak Lyman a emission 
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Figure 16(a) Calculated scattered intensity of Lyman 
a (\1216 YU) corresponding to the results shown in 
figure 12 for various values of p, with N(U)= 0.1 cm 3 
and V- 20 km sec 1 . The earth is on the axis of sym- 
metry on the upstream side of the sun (ol € = 0). For 
p< 1.0 there is a singularity at d = 180° corresponding 
to the singularities on the axis of symmetry shown in 
figures 12(a) and(b). 

occurs at 0 « 90° , and there would not be a single maxi- 
mum as observed (fig. 1 7). If p > 1 , the minimum inten- 
sity is zero according to our calculations, which are 
based on the assumption that Tj = 0. 

The observed minimum intensity of ~250 R can be 
accounted for in any of three ways: (1) the effect of 
radiation pressure is such that p < 1 ; (2) there may be a 
uniform galactic background of Lyman a [Kurt and 
Syunyaev , 1968; Tinsley, 1969, 1971; Adams ; 1971]; 
or (3) the interstellar hydrogen has a relatively high 
temperature (7/«* 10 4 °K) as discussed earlier [Thomas, 
1971 ] . If Tj is small and the effects of radiation pressure 
do not exceed those of gravity (p &O.l), then a small 
bright patch should occur near the center of the region 
of low intensity as a result of the high neutral hydrogen 
density in the vicinity of the axis of symmetry [Blum 
and Fahr , 19706] . However, the bright patch may not 



Figure 16(b) The same as (a) but with the earth at 
±90° from the axis of symmetry (Le. y oc e = 90°, 270 °) 
showing the effect of parallax. With p- 0.7, and the 
same values assumed for the various parameters involved , 
the direction of the region of maximum emission inten- 
sity can vary by about 60 ° throughout the year. 

be very prominent unless observations are made exactly 
along the axis of symmetry (a e = 6 = 180°). Preliminary 
calculations indicate that for a e ^ 180° the secondary 
intensity maximum is not at all pronounced and could 
be easily missed by observations made with wide-angle 
detectors [Johnson, 1971]. If the minimum intensity 
represents galactic Lyman a, then this should be distin- 
guishable from the scattered solar Lyman a since the line 
width is expected to be considerably greater. 
Furthermore, as pointed out by Blum and Fahr 
[1970c], galactic and local components of the diffuse 
Lyman a radiation can be distinguished by the fact that 
the former must be constant in time, while the latter 
must reflect changes in the solar Lyman a flux, espe- 
cially the 27-day variations [Meier, 1969] . Observations 
of such time variations could remove these uncertainties, 
as could observations made at greater distances from the 
sun (beyond the orbit of Jupiter) where the galactic 
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X584 Hel EMISSION RATE (Rayleighs) 



Figure 16(c) Intensity of scattered A 584 He/; we have 
used the same parameters as in (a)and(bf and assumed 
that the interstellar gas contains 8 percent He. Note that 
there is a maximum in the direction of the sun as well as 
in the direction of the downstream axis of symmetry . 
The former maximum occurs in the case of helium but 
not hydrogen because the helium is easily able to pene- 
trate to within the orbit of earth. The effect of varying 
rfEe) from 0.5 to 1.0 is shown for the case a e = 180°. 
The solar line width is assumed to be 0.023 A [e.g. t 
Donahue and Kumer, 1971] which is probably too 
small, however the scattered intensity varies approxi- 
mately inversely with the line width , and for other 
values the scattered intensity can be obtained by simply 
rescaling this diagram. 

component should be dominant. The suggestion by 
Thomas [1971] that Tj is relatively large is consistent 
with theoretical work on the intercloud component of a 
neutral interstellar gas. This is certainly a possibility if 
the background soft X-ray flux is sufficient to heat the 
interstellar gas [Silk and Werner , 1969; Werner etal. t 
1970], but probably not if it is necessary to rely on 
heating by low-energy cosmic rays. There are at least 



Figure 16(d) The same as (c ), but with a e - 270° } 
showing the displacement of region of maximum 
intensity . 



Figure 16(e) Intensity of scattered > 304 He// resulting 
from the presence of singly ionized helium in the solar 
wind for a € = 0 ° , and '‘/(He) - 0.5. Most of the scattered 
ligh t comes from the direction of the sun. The scattered 
intensity is approximately inversely proportional to the 
assumed solar line width , which in this case has been 
taken to be 0.06 A 

two ways of determining whether or not the tempera- 
ture of the interstellar gas is large: (l)from measure- 
ments of the intensity of scattered A584He/, which 
should show a very pronounced maximum in the direc- 
tion of the axis of symmetry a e = 180°, and to some 
extent for other values of a e if 7/ is small (figs. 16(c) 
and (</)); and (2) from measurements of the width of the 
scattered solar Lyman a line, made when a e = 90° , 270° , 
and away from the sun (6 - 90° , 270°), where scattering 
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should take place in the line center if Tj is small. In this 
case, a hydrogen absorption cell should have a 
pronounced effect on the observed radiation unless there 
is a substantial galactic component [. Blamont , private 
communication, 1971] . 

The resonantly scattered helium radiation (A584 He/) 
is of particular interest since helium atoms penetrate 
easily to within the orbit of earth, and hence there is 
always a maximum in the direction of the sun. The 


intensity minimum occurs in the direction away from 
the sun except when a e » 180° (figs. 16(c) and (d)). 
These effects could also be found in the case of Lyman a 
for observations made in the region r> 5 AU. In our 
calculations we have adopted a solar line width of 
0.023A [Donahue and Kumer, 1971] ; this line width 
may be too small, in which case the expected scattered 
intensity should he reduced proportionately downward 
from the values shown in figures 16(c) and ( d ). There 




Figure 17 Contour maps of the solar Lyman a background in ecliptic coordinates for 
the three spinups (SU's) of OGO 5: (a) SU 1 on 13 September 1969; (b)SU 2 on 
15 December 1969; and (c)SU 3 on 1 April 19 70 . These figures have been taken from 
Thomas and Krassa [1971]. Similar diagrams have been given by Bertaux and Blamont 
[1971]. (dj Contour map of scattered \304HeII radiation based on the observations of 
Johnson etal. [1971], and taken from Ogawa and Tohmatsu [1971]. 
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Figure 17 Continued . 


(c)SU 3 on 1 April 19 70 


have been only a few observations of the diffuse helium 
radiation to date; these indicate an upper limit of ~2R 
for the intensity of the diffuse 584A radiation from 
extraterrestrial neutral helium in a direction away from 
the sun [Young etaL,l96S\ Johnson etal., 1971; Ogawa 
and Tohmatsu, 1971]. We suggest that these observa- 
tions be pursued for various combinations of values of 
a e and 0 , and also outside the influence of the neutral 
helium envelope of the earth. 

There must be diffuse Lyman a and X304 He// radia- 
tion resulting from the presence of neutral hydrogen and 
singly ionized helium in the solar wind as indicated in 
figure 14. The Lyman a produced in this way can be 
distinguished from that described above since it should 
be red-shifted by about 1.5 A from the center of the 
solar line. However, calculations by Gregory [1971] 
indicate that for ju = 0.8, N(H) = 0.1 cm"” 3 , and a e = 0, 
the intensity of the scattered Lyman a from this source 
is only ~3.5 R in the direction 0 = 0, which is probably 
too small to be distinguished from the galactic back- 
ground. The intensity of scattered 304A radiation is 
shown in figure 16(e) for the case a e = 0; as in the case 
of radiation scattered by neutral helium, maximum 
emission occurs in the direction of the sun. For 
0 « 1 20° to 150°, the intensity of the scattered radia- 
tion is ~1 j R, which is consistent with an observation of 
1.0 ±0.2 R reported by Ogawa and Tohmatsu [1971] . 

Another possible source of diffuse Lyman a radiation 
is scattering from fast hydrogen atoms produced by 
charge exchange in the region of subsonic solar wind 
flow beyond the shock transition. In this case the emis- 
sion should be centered at ~0.3A to the red side of the 


center of the solar line, and have a half width of about 
lA as a result of the random velocities of the hydrogen 
atoms produced in this region. This emission has been 
considered by Patterson etal [1963] and Hundhausen 
[19686] in an attempt to account for an observation of 
~540 R of Doppler-broadened Lyman a reported by 
Morton and Purcell [1962], However, this observation 
showed only that there is more or less isotropic emission 
outside ±0.04A from the center of the line, which is 
easily accounted for by scattering from interstellar 
hydrogen, since it is moving at velocities exceeding 
20 km sec" 1 , and will in general produce Lyman a emis- 
sion displaced about 0.1 A from the line center. The cal- 
culations of Patterson etal [1963] and Hundhausen 
[19686] indicate that the isotropically moving hydrogen 
atoms must be produced in the region r = 10 to 20 AU 
to explain the observed emission intensity. However, 
since it seems likely that as noted, the transition to sub- 
sonic flow does not occur within 50 AU, this should not 
be an important source of diffuse Lyman a. 

It would be of interest to extend the calculations 
described in this section to include other emissions. 
Undoubtedly other hydrogen and helium lines are 
present in the diffuse radiation with intensities propor- 
tional to those of the X1216 HI and X584 He/ emissions 
[Reay and Ring, 1969; Tohmatsu, 1970] . Furthermore, 
the possibility of scattering by other species is also of 
interest, and on the basis of table 1 it would appear that 
neon in particular might produce a detectable emission, 
since r/(Ne) is small and neon is a relatively abundant 
element. 
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He 304A Isophote map projected 

into celestial coordinates (Multiply 0.36 to get 

brightness in Rayleighs) 



Right ascension 

(d) Contour map of scattered U04 Util radiation based on the observations of Johnson etal\l97l\, and taken from Ogawaand Tohmatsu [1971] 
Figure 17 Concluded . 










INTERACTION OF THE SOLAR WIND WITH 
NEUTRAL INTERSTELLAR GAS 

In this section we consider the effect on the solar wind 
of the neutral component of the interstellar gas. Various 
aspects of this interaction have been considered by 
several authors [Ax ford etaL, 1963; Patterson etal , 
1963; Ax ford and Newman , 1965; Dessler, 1967; Kern 
and Semar , 1968; Hundhausen, 1968a; Fahr f 1968 a,b, 
1969, 1970; Blum and Fahr ; 1969, \91Qa,b\ Semar f 
1970; Holzer, 1970, \91 1 \ Holzer and Ax ford, 1971; 
HM/s, 3971; Bhatnagar and Fahr, 1971] . In this discus- 
sion of the problem, we shall follow the treatment given 
by Holzer [1971] , which is similar to the work of Wallis 
[1971] but somewhat more detailed. 

It has been shown that the interstellar neutral gas 
penetrates to within a few astronomical units of the sun 
before becoming significantly affected by losses due to 
photoionization and charge exchange. Since the solar 
wind is expected to extend to a heliocentric distance of 
the order of 100 AU in the absence of any interaction 
with neutral interstellar gas, it is evident that as a first 
approximation we can treat the solar wind as if it were 
flowing through a uniform background of neutral gas 
with density N , temperature T n , and velocity V relative 
to the sun. Even with this assumption, the problem is 
very difficult since there must be nonradial flow in the 
region of supersonic solar wind flow as well as in the 
subsonic region; furthermore, the effects of the inter- 
planetary magnetic field should not be ignored. How- 
ever, a simplification can be achieved if we consider only 
the directions parallel and antiparallel to the velocity 
vector of the interstellar gas, since the flow can then be 
treated as being approximately radial. In this case, as 
noted earlier, the azimuthal component of the inter- 
planetary magnetic field is dominant in r >5 AU and 
accordingly we may assume that B r *& 0 and 
B'zB'pccfury 1 . 

If we represent the solar wind as a steady, radial, 
spherically symmetric flow of a proton-electron plasma 
through an atomic hydrogen gas, then neglecting viscos- 
ity and heat conduction, the equations of motion can be 
formulated as follows: 

i-J-(nur 2 ) = 9p (61) 


Id f 
r 2 dr L 


nurH-u 2 
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(63) 

Here c s (^\/yp/p) is the plasma sound speed, m is the 
proton mass, and c is the speed of light. The rates of 
mass, momentum, and energy production are repre- 
sented by q, 2 ^, and (f , per proton mass per unit 
volume, respectively, with the subscripts p and c refer- 
ring to the processes of - photo ionization and charge 
exchange. Thus 
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where o p is a mean photoionization cross section, & is 
the average energy of a photoelectron produced in the 
photoionization process, T(=T e + T p ) is the plasma 
temperature, otT is the proton temperature, and o c is the 
resonant charge exchange cross section. The velocity of 
the interstellar neutral gas (V) is taken to be radial. 

The electric field (E), the current (j), and the mag- 
netic field satisfy the equations 
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Thus, with the assumption that the magnetic field is 
azimuthal, we find that 


and 


7<P<BV-|fel 


Aw 


(72) 


j*E =-w 



Aw 


(73) 


We can now write the equations of motion in the form 


In the special case where B = G=V=T n = & = 0, 
equations (74) and (75) reduce to a single first-order 
differential equation in terms of the Mach number 
M = u/c s : 


M dr r r 


=3 * \ 4 4 2) 
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(74) Note that v/u is a function only of M, since in this 
special ease (69) reduces to 


^! =(r - 1) (_2 a_ (75) 
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with the following definitions: 
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^ = 0^11(1+ — - jj ) ( 81 ) 

Many features of the solutions of the more general equa- 
tions (74) and (75) can be understood from an examina- 
tion of the solutions of (80). In effect, this is the prob- 
lem considered by Semar [1970], who carried out a 
numerical integration of the equations of motion (61), 
(62), and (63), for this case taking 7=2, a - 1, and 
omitting the term enclosed in. the bracket in (81). 
Unfortunately, the existence of a singular point at M = 1 
was overlooked in Semar’ s analysis, which consequently 
implies that a shock-free (“transonic”) transition to sub- 
sonic flow will always occur. In fact, as pointed out by 
Wallis [1971] and Holzer [1971] , this is not necessarily 
the case. 

Equation (80) has a singular point at which (M 2 - 1) 
and the right-hand side of (80) vanish simultaneously. If 
this point is defined to be at r = r c , then 

^ = ^{3-7 -2Rf c 

+ [ AR 2 f c 2 + (20 + Ay)Rf c + (y-3) 2 ] 1 /2 } (82) 
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where, using (64) and (69), and assuming that the 
neutral gas is optically thin in the region r Q <r<r c , we 
have 


where the subscript o denotes quantities evaluated at a 
point r~r Q close to the sun. 
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Note that under typical solar wind conditions and with 
0.05 < N < 1.0 cm" 3 , r c « 1/A. 

The form of the characteristics passing through the 
singular point can be determined by expanding (80) 
about M = 1 , r - r c . Defining X = M - 1 and £ = (r/r c )-l , 
we obtain to first order in X and £ 



■RADIAL DISTANCE (r/r c ) 


(a) b 2 + 4a >0 



(88) (b)b 2 +4a<0 


A characteristic passing through the singular point has a 
slope at the singular point given by S = lim(X/£). Hence S 

I 

can take one one of the two values 

8+ = ~ (b ± \/b 2 + 4a ^ (89) 

If b 2 + 4a > 0] the characteristics passing through the 
singular pointtcan have slopes S + or S_. If b 2 + 4a = 0, all 
characteristics passing through the singular point have 
slope S + = S_~ 6/2. If b 2 + 4a < 0, the characteristics in 
the vicinity of the singular point are spirals and the 
corresponding solutions are multivalued in M. A family 
of solutions of (80) ■ for which b 2 + 4a > 0 is shown in 
figure 18(a) and a family for which 6 2 + 4a < 0 is shown 
in figure 18(6). 


Figure 18 Characteristics of equation (58) [from 
Holier, 1971 ] . Note that a transonic solution is possible 
in (a) but not in (b) where the singular point is a focus. 

On examination of the possible solutions displayed in 
figures 18(a) and (6) one finds that only one character- 
istic leads to infinity in each case, and that these are 
such that the pressure p **> 0 as r <*>; there are no solu- 
tions that yield a finite pressure at infinity . For the case 
of moderate and low interstellar densities, the character- 
istics are as shown in figure 18(a), and it is always 
possible to find a solution starting with suitable condi- 
tions at r = 5 AU, say, which contains a smooth transi- 
tion to subsonic flow and eventually has vanishing 
pressure at infinity. For somewhat higher interstellar 
densities (N> 10 cm" 3 , depending on the value chosen 
for other parameters), the characteristics are likely to be 
as shown in figure 18(6); in this case, again it is possible 


638 





to find a solution starting with suitable conditions at 
r = 5 AU, and with vanishing pressure at infinity, but it is 
necessary to insert a shock at a correctly determined 
location near the singular point to make the solution 
single valued. In the numerical treatment carried out by 
Semar [1970] these complexities have passed unnoticed, 
although the results obtained are probably a fairly good 
representation of the solution except in the vicinity of 
the singular point. It should be noted, however, that the 
form for the collision frequency v used by Semar is 
inadequate if the Mach number is not large, and the 
former has an important effect on the properties of the 
characteristics in the subsonic region. 

The solutions of equation (80) are useful in that they 
suggest what we can expect in the more general case 
involving a pair of differential equations — namely, (74) 
and (75) — that must be solved simultaneously for two 
unknowns u andc s . Holzer [1971] has treated this prob- 
lem by a straightforward integration of (74) and (75), 
having assigned values for u , p, T, and B at a point 
relatively close to the sun (r = 5 AU) where the solar 
wind can be expected to be essentially unaffected by the 
presence of neutral interstellar gas. The integration 
proceeds in the direction of increasing r, until the Mach 
number Mj- uj(c s 2 + c^ 2 ) 1 ' 2 decreases to unity, or (if 
a shock transition is inserted) until the solutions 
approach some asymptotic form that can be recognized. 
In this procedure, there is some danger that singular 
points might be missed and that the solutions might be 
unstable; however, the results obtained appear to be 
satisfactory. Examples of solutions of equations (74) 
and (75) obtained in this manner are shown in fig- 
ures 19(a) and ( b ) and 20(a)-(c). The initial conditions 
have been taken to be u 0 = 400 km sec” 1 , 
n Q = 0.2 cm" 3 , c so = 15 km sec” 1 , B 0 = 0.7X10“ 5 gauss 
at r = r Q = 5 AU, together with a range of values for the 
density (TV) and velocity (V) of the interstellar neutral 
gas. 

The location of the sonic point (M ~ 1; r = r c ) 
is shown in figure 19(a) for a range of values of AT and V; 
this yields an outer limit for the position of any possible 
shock transition in the flow. Note that the location of 
the sonic point is quite sensitive to the density of the 
neutral interstellar gas, varying roughly as AT 1 . Further- 
more, there is a substantial difference between the values 
of r c on the upwind (F<0) and downwind (K>0) 
sides of the sun, which should produce a distinct asym- 
metry in the region of supersonic solar wind flow — for 
example, for N~ 0.1 cm" 3 and V- 20, -20 km sec” 1 , 
r c = 250 AU, 400 AU, respectively. The variation of 
solar wind speed with distance is shown in figure 19 (b) 
for K~-20kmsec" 1 and various values of N. If we 



(a) Radial distance to the sonic point (see 
figs. 16(a), (b)) for various assumed values of the relative 
velocity of the interstellar gas , and interstellar hydrogen 
densities in the range 0.0 to 1.0 cm 3 . 



(b) Solar wind velocity as a function of radial distance 
from the sun for the case in which the interstellar gas has 
a velocity of 20 km sec 1 towards the sun , and for var- 
ious assumed values of the interstellar hydrogen density. 
Note that for the preferred value (0.1 cm 3 ) the solar 
wind velocity is reduced to < 300 km sec 1 at a distance 
of 100 AU where the shock transition might be expected 
to occur (see fig. 7). 

Figure 19 Effect of interstellar neutral gas hydrogen 
on the solar wind [ from Holzer, 1971] . 

assume that N= 0.1 cm" 3 and V = -20 km sec" 1 , as 
suggested by the observations described in the preceding 
section, it can be seen that the solar wind speed is 
reduced from 400 km sec” 1 at 5 AU, to about 
270 km sec" 1 at 100 AU, which we have estimated to be 
the likely position of the shock transition. Since photo- 
ionization is a relatively unimportant process the mass 
flux in the solar wind is essentially unaffected by the 
deceleration and hence the momentum flux (or ram 
pressure) is reduced in proportion to the velocity. This 
result can also be deduced easily by noting that the 
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MACH NUMBER (M^u/fc^cjJ)''*) 



(a) Variation of the Mach number with radial distance 
for various assumed positions of the shock transition . 
Note that the solutions terminate at a finite distance 
from the sun if the shock is inserted too close to the 
critical poin t 



(b) Variation of the plasma density n, total pres- 
sure pji , and magnetic field pressure B 2 j 'Sir, as a 
function of radial distance for the case where the shock 
is inserted at 15 A U. 



(c) Variation of the solar wind speed u, the sound speed 
c s , and the Alfveh speed cj, with radial distance for the 
case in which the shock transition is inserted at 15 AU. 


Figure 20 Variation of solar wind parameters with distance for the case 
V= 20 km/sec ~ l , and N= LOcmT 3 [from Holier, 1971]. Note that the solutions are 
somewhat similar to those shown in figure 7 for distances less than a few hundred AU. At 
very large radial distances the solutions may not be at all realistic; however , it is inter- 
esting that the asymptotic forms create a balance between the stresses associated with the 
magnetic field and the frictional interaction between the plasma and the neutral inter- 
stellar Wind. 
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mean free path of solar wind protons against charge 
exchange is approximately 300 AU if N= 0.1 em“ 3 and 
hence on reaching a heliocentric distance of 100 AU, the 
plasma must lose (l-e“ 1A3 ) « 1/3 of its momentum per 
unit mass. Formally, this is equivalent to making the 
approximation Mj 2 » 1 in (74) and putting 
=“[( 7 + \)j2\w. Thus we see that loss of 
momentum due to charge exchange within the region of 
supersonic solar wind is likely to have a significant effect 
on the location of the shock transition (in the sense that 
r s might be reduced by 10 to 15 percent from the value 
estimated previously. 

Solutions of equations (74) and (75) for the 
postshock flow have also been considered by Holzer 
[1971], again with the assumption that the flow is 
strictly radial; an example with F = 20 km sec ” 1 and 
N= 1 cm ” 3 is shown in figures 20(a)— (c). There appear 
to be no solutions that extend to infinity in the upwind 
direction (V < 0); in the downwind direction (F> 0), 
however, such solutions exist, provided the transition to 
subsonic flow takes place within a certain critical radius 
(r=rj). The Mach number is shown as a function of 
radial distance in figure 20(a). Note that if the shock 
location is such that r s >r x , the solution in the subsonic 
region bends upward, intersecting the line Mj- 1 with 
infinite slope, and does not extend to r^°°. In cases 
where r s <r t the solutions in the subsonic region 
asymptotically approach the same curve, which is of the 
form Mjr 312 = constant. The behavior of the other 
dependent variables for the case = 15 AU <r l is 
shown in figures 20(b) and (c). The flow in the subsonic 
region can be divided into three parts: ( 1 ) in 
15 AU 0< 50 AU, the plasma is hot but is rapidly 
cooled as a result of charge exchange, thus causing the 

velocity to decrease and the magnetic field strength to 
increase; (2) in 50 AU < r< 10 3 AU, the magnetic field 
dominates the total pressure, but the velocity is con- 
trolled by photoionization and charge exchange and is 
approximately equal to F; (3) in r> lO 3 AU, the mag- 
netic field strength increases and the velocity decreases 
so that the asymptotic state is controlled by a balance 
between the tension stresses of the magnetic field and 
the frictional drag due to photo ionization and charge 
exchange. 

It is evident from the scales involved that we do not 
have to take these last results too seriously. Nevertheless, 
they do demonstrate rather clearly the complexity of 
the interaction that can be expected in more realistic 
situations when nonradial flow is important. We see, for 
example, that in the subsonic region on the upwind side 
of the heliosphere the most important effect of charge 


exchange is to cool the plasma, thus permitting the mag- 
netic field strength to grow more rapidly than it would 
otherwise (fig. 8 ). Charge exchange should not be impor- 
tant in the upwind subsonic region as far as the 
pressure/ momentum balance is concerned since the 
streaming velocities are low. However, the increase of 
magnetic field strength enhances the inward force due to 
the tension stress associated with the field, and as 
noted earlier, this could produce a small but significant 
radial pressure gradient in the subsonic region. This 
effect, together with the loss of solar wind momentum 
due to charge exchange, could easily render our previous 
estimate of r s = 100 AU for the location of the shock 
transition too large by perhaps 20 to 30 percent [Kern 
and Semar, 1968] . However, a better self-consistent 
treatment is required before we can be sure that this is 
the case. 

On the downwind side of the sun (in the “tail” of the 
heliosphere) we can expect the solar wind in the super- 
sonic region to lose momentum due to charge exchange, 
so that the shock transition (if any) is not as far from 
the sun as it would be otherwise. In the subsonic region, 
the plasma should at first be cooled by charge exchange 
until drops to such low values that the magnetic field 
has an important influence in controlling the flow. 
Momentum exchange due to photoionization and charge 
exchange will maintain the velocity of the flow at 
approximately that of the interstellar neutral gas. Simi- 
larly, the temperature of the plasma should ultimately 
be controlled by the input from photoelectrons and by 
the temperature of the interstellar gas. There is no strong 
reason for believing that the magnetic field will eventu- 
ally increase as indicated in figure 20(b) since at the 
radial distances involved (r> 10 3 AU) in the actual 
heliosphere the magnetic field will certainly not be 
azimuthal as assumed in the calculations. Indeed, it 
should be expected that at distances of the order of a 
few hundred astronomical units the tail of the helio- 
sphere should gradually deflate and collapse as the 
thermal pressure is removed from the plasma as a result 
of charge exchange. As the (3 of the plasma decreases to 
small values the magnetic structure of the tail should 
disintegrate due to field line reconnection at any sur- 
viving neutral sheets. Eventually, the field must become 
completely connected to the interstellar magnetic field 
and the plasma will disperse into the interstellar medium 
along the field lines. Thus we do not expect the tail of 
the heliosphere to be very long; indeed, after the first 
few hundred astronomical units, the term “wake” would 
be more appropriate. Finally, it should be noted that as 
a result of the very low densities involved, the plasma in 
the tail and wake is unlikely to be rapidly neutralized by 
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recombination unless the temperature of the interstellar 
neutral gas is relatively low (Tj *** 10 2 °K). 

THE INTERACTION BETWEEN GALACTIC 
COSMIC RAYS AND THE SOLAR WIND 

Although galactic cosmic rays must pass through the 
outer regions of the heliosphere to reach the earth, they 
do not appear to carry much information about these 
regions. Nevertheless, it is useful to consider various 
aspects of the behavior of galactic cosmic rays within the 
heliosphere, especially as it is likely that the situation 
will change when observations from deep space probes 
become available. 

At very high particle energies (> 200 GeV), it is 
found that the solar wind and the interplanetary mag- 
netic field have no effect on the mean cosmic ray inten- 
sity observed at the earth in the sense that there is no 
solar cycle modulation or evidence for Forbush 
decreases. It might be expected, therefore, that if the 
cosmic ray “gas,” including these particles, has a velocity 
V C r relative to the solar system, then this should be 
detectable as a sidereal diurnal variation of the intensity. 
That is, the differential intensity in the kinetic energy 
range (T, T+ dT) should have the form 

KT,t)=] 0 (mH(T)cosm] (90) 

where t is the time, <p(t) is the angle between the 
detector acceptance direction and and %{T) is the 
anisotropy given by 


£(D = 


3CV cr 


(91) 


where v is the particle speed corresponding to kinetic 
energy T , and C is the Compton-Getting factor: 


aX)=l ~h w { olTU) (92) 


where a = (T+ 2E 0 )I(T+ E Q ), E Q is the particle rest 
energy, and U( T) = 47r/(7)/v is the differential number 
density in (T, T + dT) [Gleeson and Axford, 1968# ; 
Forman , 1970] . At high energies (T» E Q ), it is a good 
approximation to take a =1 , v=c , and ]{T) 
where ju - 2.5. Thus C— (2 + ju)/3 — 1,5, and if 
V cr “ V — 20 km sec“ 1 , we should expect that 
| — 0.03 percent, in principle an observable value. 

Unfortunately, the situation is not as simple as this 
analysis suggests, since although the irregular component 


of the interplanetary magnetic field does not signifi- 
cantly affect particles with energies of the order of 
100 GeV and higher, the regular (Archimedes spiral) 
field can bend the particle trajectories so that any anisot- 
ropy becomes obscured. This is evident in figure (21), 
which also shows contours on which particles have gyro- 
radii in the local mean interplanetary magnetic field 
equal to the scale length of the field (B/VB). Within 
each contour, the trajectories of particles with energies 
less than the designated energy can be completely turned 
by the magnetic field. McCracken [unpublished, but see 



Figure 21 Contour diagrams similar to those shown in 
figure 3 indicating regions of the interplanetary medium 
in which cosmic rays of a given energy have a gyroradius 
in the local magnetic field which is less than the charac- 
teristic scale of the field (B/VB). The diagram is some- 
what misleading in the low latitude region beyond the 
few astronomical units since it does not take into 
account the sector structure of the interplanetary mag- 
netic field. Since the sectors have a characteristic scale of 
only a few astronomical units , the trajectories of parti- 
cles with gyroradii much larger than this will be 
essentially unaffected by the field. 


Antonucci et al, 1970] has examined the nature of the 
sidereal variations to be expected in the energy range 
50—500 GeV (the range appropriate to underground 
muon detectors). He has shown that diurnal and semi- 
diurnal intensity variations are to be expected, with the 
diurnal variation varying in amplitude throughout the 
year in such a manner that sidebands in the modulation 
should occur at 364 and 368 cycles/yr. 
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There is some evidence for the existence of a sidereal 
anisotropy obtained from detectors at 60 to 80 m.w.e. 
[Elliot et al, 1970; Antonucci et al, 1970] ; however, it 
is not wholly convincing. For example, underground 
muon detectors at depths less than about 165 m.w.e. are 
sensitive to particles having energies less than 100 GeV 
[Ahluwalia, 1971], and the solar diurnal modulation 
(to which such particles are clearly subject) can produce 
pseudosidereal effects that can contaminate any genuine 
sidereal variation [Swinson, 1971] . It is probably neces- 
sary to make use of extensive air shower observations to 
obtain a good value for the cosmic ray anisotropy, since 
the particles involved have energies in excess of 10 12 eV 
and should not be noticeably affected by the inter- 
planetary magnetic field or solar modulation. Unfortu- 
nately, the intensity of cosmic rays with energies 
> 10 1 2 eV is low, and so far it has been possible only to 
place an upper limit of about 0.1 percent on the 
anisotropy. 

At energies less than 100 GeV, galactic cosmic rays 
are modulated in intensity in a variety of ways by their 
interaction with the solar wind and interplanetary mag- 
netic field \ Jokipii [1971] has reviewed the theory of 
some of these effects. In particular, the cosmic rays 
undergo a solar cycle modulation that appears to be the 
result of a reduction of intensity in the inner solar sys- 
tem due to outward transport of the particles by the 
solar wind with associated energy changes [Parker, 1963 , 
19656; Gleeson and Axford , 1967; Gleeson , 1969; 
Jokipii and Parker , 1970] . The equations describing the 
behavior of the differential number density U and 
current density S, assuming that the solar wind is steady 
and radial, and taking only radial diffusion into account, 
are 


Urch [1971] have been helpful; however, for many 
purposes it is more convenient to use various approxi- 
mate forms of equations (93) and (94) such as those of 
Fisk and Axford [1969] and Gleeson et al, [1971] that 
yield asymptotically valid solutions. Gleeson and Axford 
[1967, 19686] and Fisk and Ax ford [1969] have shown 
that an adequate approximation for T ^ 200 
Me V/nucleon at the earth can be achieved by neglecting 
the term S on the left of equation (94); the resulting 
equation 


CuU=k ^ (95) 

dr 

is termed the force-field equation. 

If k(t, T) is a separable function of r and T — that is, 
k(t, T) = Ki(r)K 2 fT) - this equation can be integrated to 
yield a result similar to Liouville’s theorem: 

j(r,E) _ j(°°,E + $) ^ 

E2 ~ E o 2 (E + <f>) 2 -E 0 2 

where E~T+E 0 is the total energy of the particles, and 
d> can be regarded as a “potential” associated with the 
(apparent) force exerted on the cosmic ray distribution. 
In particular, for the case in which k 2 a Pv/c = al V/c, 
which appears to be a reasonable choice in the vicinity 
of the earth for rigidities P>lGv [Jokipii and 
Coleman , 1968] , = lZel0(r), where 


#0 = 



u(r) 


dr 


(97) 


££***>•-£. j&WW (93) 

a u 

S = CuU-k — - (94) 

or 

where u is the solar wind speed and k(t, T) the effective 
radial diffusion coefficient for cosmic rays in the 
(irregular) interplanetary magnetic field. 

It is rather difficult to obtain solutions of equa- 
tions (93) and (94) with the boundary conditions appro- 
priate to the modulation problem ( U -» £/<*,( 7) as r -> 
f i S -*Q as r-* 0), although classes of exact solutions 
exist if it is assumed that a can be treated as a constant 
[Fisk and Axford , 1969] . Purely numerical treatments 
of the problem such as those of Fisk [1971a] , Gleeson 
and Urch [197 16] , Lezniak and Webber, [1971], and 


It is found that equations (96) and (97) provide a very 
good description of the (11-yr) solar modulation of 
protons, a particles, and heavier nuclei for energies 
greater than 100-200 MeV/nucleon. The “potential 
difference” between the earth and infinity is typically 
found to be in the range 150—500 MV, and the modula- 
tion shows the predicted Z dependence [Gleeson and 
Axford, 19686; Lezniak and Webber, 1971]. Goldstein 
et al [19706] and Gleeson and Urch [1971a] have 
found from numerical solutions of the full equa- 
tions (93) and (94) that galactic cosmic rays with initial 
energies less than 100-200 MeV/nucleon are effectively 
excluded from the interplanetary medium at the earth’s 
orbit (fig. 22), which is consistent with the idea that the 
modulation can be described in terms of a force field. 
Clearly it should be considered a prime task of any deep 
space mission to measure the spectra of various species 
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Figure 22 A series of essentially monoenergetic proton 

spectra in interstellar space (solid line) and the corre- 
sponding modulated spectra at earth (dashed lines) 
obtained from numerical solutions of equations (93) and 
(94) [from Goldstein etal , 1970b], The “force-field” 
solution (eq. (96)) provides a very good fit to the 
envelope of the modulated spectra for energies above 
about 10 2 MeV It is important to note that particles 
with energies <200 MeV in interstellar space (i.e., 
unmodulated peaks 1-3) make a negligible contribution 
to the spectrum at earth , where this energy range is 
dominated by particles having higher energies in inter- 
stellar space (i.e., unmodulated peaks 4-6) and which 
have lost energy as a result of the modulating process. In 
effect , galactic cosmic rays with energies less than 
100 MeV are unobservable at earth. 

of cosmic rays as a function of radial distance from the 
sun as a basis for determining the unmodulated spectra 
and the variation of the diffusion coefficient with 
distance. 

The radial gradient of the cosmic ray density or inten- 
sity at energies greater than ~200 MeV/ nucleon is given 


directly by (97) as 

1 &U I 3/ _ Cu 

U dr ~ j dr k (98) 

On taking a value for the diffusion coefficient consistent 
with that obtained by Jokipii and Coleman [1968] , one 
finds that the gradient is of the order of 
20-30 percent/ AU for 1 GeV protons, and decreases at 
somewhat higher energies approximately as 1 /P. There is 
some evidence that gradients of this general magnitude 
exist, at least for P> 1 GV, from spacecraft observations 
[O'Gallagher and Simpson , 1967; OGallagher, 1967], 
measurements of radio activity in meteorites [Forman et 
al , 1971] , and observations of the associated cosmic ray 
anisotropy perpendicular to the ecliptic [Yoshida etal., 
1971]. 

If such gradients were to exist throughout the region 
in which the solar wind is supersonic (say, r < 50 AU) 
then the total energy density of galactic cosmic rays 
outside the heliosphere would be enormously larger than 
that observed at the earth. In fact, it seems unlikely that 
the unmodulated spectrum for T> 1 GeV has a steeper 
slope than a spectrum with jffyccT^ with ju — 2.5. 
This spectrum matches the (apparently unmodulated) 
observed spectrum in T>10 2 GeV. It would require 
that the energy density of galactic cosmic rays be greater 
by only a factor of 2 or 3 outside the heliosphere than 
near the earth for T> 1 GeV, arid hence that the radial 
gradient cannot be maintained at the value quoted above 
for more than a few AU from the sun. This argument is 
essentially the same as that of Gleeson and Axford 
[19686], who for convenience assumed that 
Ki(r)cc / lr o i n anc | f oun d that if the “potential 
difference” between the earth and infinity is 
100-200 MV, then r Q must be of the order of 1 AU. 

Additional confirmation that the modulation of 
galactic cosmic rays with T> 1 GeV occurs in a region 
that is effectively contained within a few AU of the sun 
arises from comparisons of the solar cycle variation of 
the cosmic ray intensity (measured by neutron monitors, 
for example) and various indices of solar activity 
[Simpson, 1962; Simpson and Wang, 1967, 1970; 
Hatton etal. , 1968; Guschina etal., 1970; Kolomeets 
etal, 1970; Wang, 1970]. It is argued that the lag 
obtained from a cross correlation analysis involving 
these quantities should provide a rough estimate of the 
time required for the solar wind to traverse the modu- 
lating region. A cross correlation between cosmic ray 
intensity and the geomagnetic activity index Kp yields a 
lag of the order of 6 to 12 months, which suggests that 
the radius of the modulating region might be as large as 
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50 to 80 AU [Dorman and Dorman , 1966] . However, it 
has been found that if a cross correlation is made 
between the coronal green line (X Fe XIV 5303), the lag 
is only of the order of 1 month, which yields a radius of 
~7 AU for the modulating region of particles detected 
by neutron monitors [Simpson and Wang, 1967 ; Barker 
and Hatton, \91Q\Kolomeets etal, 1970] . 

It is not possible to draw similar conclusions for the 
low energy end of the galactic cosmic ray spectrum, 
simply because the degree of modulation at low energies 
is expected to be so large that we are completely igno- 
rant of the form of the unmodulated spectrum (fig. 21). 
Furthermore, at the lowest energies, the spectrum 
observed near the earth appears to consist mostly of 
particles of solar origin [Ax ford, 1970] . The effects of 
adiabatic deceleration associated with the expansion of 
the solar wind dominate the behavior of these particles. 

It seems likely that galactic cosmic ray protons and 
heavier nuclei with energies of the order of 
1—10 MeV/nucleon are to a large extent excluded from 
the region of supersonic solar wind flow, although they 
may penetrate the region of subsonic flow more easily 
since the effects of convection in this region are less 
important. At the shock transition separating these 
regions it is possible for significant acceleration of low- 
energy particles to occur, thus producing, in a sense, a 
local ‘‘source” of cosmic rays [Jokipii, 1968]. This 
effect can be demonstrated by means of appropriate 
solutions of equations (93) and (94). Let us suppose that 
the effective radial diffusion coefficient is small 
(k « ur s ) and that the solar wind convects energetic 
particles of solar origin toward the shock in such a man- 
ner that U(r,T) U 0 (T) = for (r s -r)u/a In 

these circumstances, the effects of spherical divergence 
can be neglected in r<r s , and we find that equa- 
tions (93) and (94) can be replaced by 

S = u 1 U 0 -l Ul -^:(aTU) (99) 


S = ui |jj 

expj^. (r-r s )J J T* (102) 

where B is a constant obtained in the integration of 
(100), and it is assumed that U TV everywhere. 
Beyond the shock transition (in r> r s ), we can treat the 
flow as if it were incompressible to a first approxima- 
tion, so that u = (l/4)w 1 (r s /r) 2 ; hence, equations (93) 
and (94) become 

<103) 

luAU-C(T)](jj % ( 104 ) 

where C(T)= CT ~ ^ arises from the integration of (93). 
The appropriate solutions of these equations are such 
that U, S^Oasr-^ 00 : 



Provided there is no source of particles within the shock 
itself, both U and S must be continuous at the shock 
[Gleeson and Ax ford, 1967] , and hence the constants B 
and C can be determined: 

B , „ M[l-W3Xy-l)l c - 3 “' l > 4 (107) 

oc(jx- 1) 4K 2 a(ji-\) 


— (loo) 

having put u = u l and (3/3 r)(r 2 S)^r 2 (dS/dr), and inte- 
grating (93). Thus we find that if k = tc x is assumed to be 
independent of r and T , and a is taken to be constant, 
then 

U= |^4 + Z?exp£— | (101) 


Note that the enhancement of the cosmic ray density 
resulting from the presence of the shock (C/A) depends 
on the value of the diffusion coefficient in the subsonic 
region, and may be very large if K 2 /u l r s « 1. The 
cosmic ray intensity near the earth is not significantly 
affected by the presence of this enhancement since if 
K«ur the particles are convection dominated and 
downstream boundary conditions have little influence. It 
is possible to solve equations (93) and (94) exactly with 
the solar wind distribution assumed in the above analy- 
sis, provided a is taken to be constant, k is independent 
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of T, and U<^T~^ [Fisk, 1969] . The solutions are con- 
sistent with those obtained here in the limit ur s /K » 1 ; 
they also show that it is possible to produce an enhance- 
ment of the intensity of galactic cosmic rays in the 
vicinity of the shock transition under suitable condi- 
tions. On this basis, therefore, it should be expected that 
a low-energy cosmic ray enhancement similar to those 
associated with propagating shocks [Fisk, 1971b] is 
associated with the shock transition terminating the 
supersonic solar wind, and this may provide useful 
confirmation of the presence of the transition. 

The possibility that the solar wind itself might be 
significantly affected by cosmic rays has been considered 
by several authors [Ax ford, 1965 ; Ax ford and Newman, 
1965; Dorman and Dorman, 1968; Modisette and 
Snyder, 1968; Sousk and Lenchek, 1969; Wallis, 1971]. 
In the early analyses it was assumed that there is no 
exchange of energy between the solar wind and the 
cosmic rays, and that the cosmic ray pressure 7r satisfies 
the convection-diffusion equation 

di t 

K — - un 

dr (108) 

where K(r) is the “effective” diffusion coefficient and 

it = - r otTUdT (109) 

3 J o 


In more recent treatments [Sousk and Lenchek, 1969; 
Wallis, 1971} , which take the energy exchange into 
account, equation (108) is used together with the 
following equations For the conservation of mass, 
momentum and energy of the solar wind plasma: 


d_ 

dr 


(pur 2 ) * 0 


( 110 ) 


du _ dp GTtjp dn 
pU T~~!f- r 2 ~dk 


( 111 ) 


1 d 
P" dr 


\pur* 



£ 

9 



«£i an) 

dr 


where / °° ST dT is the net radial flux of energy in 
o 

cosmic rays; since udn/dr > 0, it is evident that there is a 
net transfer of energy from the solar wind to the cosmic 
rays [cf. Jokipii and Parker, 1967] . 

In the treatment of the problem by Axford and 
Newman [1965] and others, the right-hand side of (1 12) 
was omitted; since this term represents an effective heat 
sink for the solar wind, the results obtained by these 
authors are such that the Mach number decreases toward 
unity too rapidly with increasing radial distance. Wallis 
[1971] has pointed out that when equa- 
tions (109)— (112) are reduced to a single equation for 
the Mach number in terms of radial distance, the 
characteristics of the equation have properties similar to 
those described in the preceding section for the case of 
interaction with interstellar neutral gas. In particular, 
there is a singular point at a large radial distance, the 
nature of which depends on the values assumed for the 
various parameters involved. It is found, however, that if 
the interaction is sufficiently strong that a significant 
effect is produced on the solar wind, then the cosmic ray 
pressure required in the interstellar medium is unaccept- 
ably high (more than 10 times the value observed at the 
position of the earth). An example of a solution 
obtained by Sousk and Lenchek [1969] for the case 
where K is independent of r is shown in figure 23; it is 
evident that quite unreasonable cosmic ray pressures are 
necessary in this case to reduce the solar wind speed by 
as much as 25 percent. As a corollary to calculations of 
this sort, it is found that the cosmic ray pressure cannot 
be prevented from becoming excessive unless the net 
diffusion coefficient K increases more strongly than 
linear with radial distance. This is consistent with the 
conclusion of Simpson and Wang [1967] and others that 
the dimensions of the modulating region must be rela- 
tively small (~7 AU). 

Cosmic rays are affected by the regular component of 
the interplanetary magnetic field as well as by irregular- 
ities, and as a consequence they tend to rotate with the 
magnetic field about the sun. The resulting anisotropy is 
observed at the earth as a solar diurnal variation of the 
cosmic ray intensity, which is of the order of Q.3 to 
0.4 percent, given approximately by 


On combining equations (1 10)— (1 12), we can show that 
/?/p7 - constant. The right-hand side of equation (1 12) is 
easily shown from (93) to be given by 



,22iSi! aw) 

V 

where £2 is the angular velocity of the sun [Ahluwalia 
and Dessler, 1962; Parker, 1964; Axford, 1965]. This 
anisotropy has yet to be detected at low energies 
(T<2 GeV) from spacecraft; it is possible that it is 
absent at energies less than about 200 MeV/nucleon 
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Figure 23 Variation of the solar wind speed V, sound 
speed Q and cosmic ray pressure n, with distance from 
the sun for the case in which the effective diffusion 
coefficient K is constant [from Sousk and Lenchek, 
1969] . It should he noted that in order for the solar 
wind speed to drop to 300 km se<T x , the cosmic ray 
pressure must rise by a factor of 10 over its value at the 
orbit of earth, which seems unacceptably large . The solu- 
tions have not been carried through the singular point of 
the equations where the solar wind and the sound speed 
are equal [cf. Wallis , 1971 ] . 


where C(7)^0 [Gleeson and Ax ford, 1968a]. It is 
difficult to escape the conclusion that this anisotropy 
must increase with increasing heliocentric distance, and 
become very large indeed at distances of the order of 
50 AU. However, the anisotropy can be suppressed in 
two ways: by (1) a gradient of the intensity perpendic- 
ular to the ecliptic, and (2) isotropy of the diffusion 
tensor. The enormous intensity gradients necessary for 
any significant effect in this respect are difficult to 
accept, but the second possibility may well be impor- 
tant. A more accurate form of equation (1 14) is 


fyC T) = 


3 C£lr sin 0 


K||- 


“1 


k || + k tan 2 i 


(115) 


where k(| and /cj_ are the diffusion coefficients parallel 
and perpendicular to the mean magnetic field direction, 
respectively. It has been pointed out by Jokipii [19 66b, 
1971] and Jokipii and Parker [1969] that as a result of 
the wandering of magnetic field lines with respect to the 
direction of the mean field it is possible for in 


which case the anisotropy given by (115) must dis- 
appear. This effect is probably not as pronounced at the 
orbit of earth as suggested by Jokipii and Parker [1969] 
[Subramanian, 1971] ; however, it is expected to 
increase in importance with increasing heliocentric 
distance and ultimately may be capable of keeping ^ 
within reasonable limits. A breakup of the interplanetary 
magnetic field sector structure as described earlier would 
serve the same purpose. It should be noted that in the 
presence of a large value of the medium comprising 
the solar wind plasma and the cosmic ray gas would be 
liable to streaming instabilities [Wentzel, 1969], which 
in turn would tend to reduce both k || and ki and thus 
increase the overall modulation. 

The problem of how cosmic rays penetrate the region 
occupied by the solar wind and interplanetary magnetic 
field is of some interest. In the absence of direct connec- 
tion between the interplanetary and interstellar magnetic 
fields it is difficult to understand how low-energy elec- 
trons, for example, could enter the interplanetary 
region with any reasonable efficiency. Parker [1968] has 
suggested that the effect of field line wandering 
described above might be adequate to bring any given 
field line sufficiently close to the interface between the 
solar plasma and interstellar region to allow it to receive 
its full quota of cosmic rays, regardless of where the 
field line intersects the surface of the sun. Alternatively, 
we would argue that the access of cosmic rays can take 
place most easily if the interplanetary field is connected 
directly to the interstellar medium. There is little reason 
for believing that this is not so, and in the case of the 
interplanetary magnetic field and the geomagnetic field, 
for example, the evidence for interconnection is very 
compelling [Akasofu and Ax ford, 1971] . 

Schatten and Wilcox [1969, 1970] have suggested 
that interconnection of the interplanetary and inter- 
stellar magnetic fields could explain the 20-year cycle of 
the diurnal variation of the cosmic ray intensity at the 
earth reported by Forbush [1967, 1969]. It is argued 
that the effect is related to the 20-year solar magnetic 
cycle, and that the access of cosmic rays is controlled to 
some extent according to whether or not the solar mag- 
netic field is favorably directed so that cosmic rays have 
easy access to the heliosphere. In fact, the situation is 
quite complex, and it is not obvious that this is the 
correct explanation. However, it should be noted that a 
similar situation exists in the case of access of aniso- 
tropic fluxes of solar energetic particles into the earth’s 
magnetosphere, where the direction of the inter- 
planetary magnetic field can lead to a pronounced 
asymmetry in the particle fluxes observed at low alti- 
tudes in the polar magnetosphere [Reid and Sauer, 
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1967; Engelman etal, 1971; Van Allen etal , 1971]. 
The possibility that interconnection of the interplane- 
tary and interstellar magnetic fields might affect the 
intensity of low-energy cosmic ray electrons has been 
discussed by Fisk and VanHollebeke [1971] in an 
attempt to explain the “quiettime increases” of these 
particles [Simnett etal , 1971]. Kovar and Dessler 
[1967] have suggested that the sidereal anisotropy of 
galactic cosmic rays might be related to the asymmetry 
of the outer heliosphere. As noted previously, however, 
the particles that can be expected to show a sidereal 
anisotropy must have such high energies that they pass 
unaffected through the complex field structure of the 
outer heliosphere. As a result of the configuration of the 
interplanetary magnetic field, any effect associated with 
the asymmetry of the heliosphere will lead to spatial 
variations of the mean intensity of cosmic rays rather 
than a sidereal anisotropy. 

APPENDIX: A HYDRAULIC ANALOGY TO THE 
SOLAR WIND 

Consider the flow resulting from a jet of liquid incident 
on a horizontally held dinner plate. Under steady 
conditions, the total mass flux 

pQ = litprhu (Al) 

is constant, where p is the density of the liquid, r is the 
radial distance from the center of the plate, h(r) is the 
depth, and u(r) the mean radial velocity of the liquid. It 
can be shown that the flow pattern divides in general 
into a supercritical region (r < rj) and a subcritical region 
(rj<r<R = radius of the plate), in which the Froude 
number F = u/(gh) l/2 is greater than or less than unity, 
respectively. 

Conservation of radial momentum requires that 

2 + i gh ^ h = constant = ~ gH 2 (A2) 

where H is the “total head.” Accordingly, if F» 1 , u is 
approximately constant and h(r) \jr\ if F« 1, h is 
approximately constant and u a 1/r. In the case of the 
solar wind, where the Mach number plays the same role 
as the Froude number, these results are analogous to the 
constancy of u when M » 1 , and of p when M « 1 . 

The regions of supercritical and subcritical flow are 
separated by a hydraulic jump at r-.fy, across which 
conservation of mass and momentum require that 

Mi =h 2 u 2 = ~ (A3) 

2717 

^! 2 +lgh hi = (u 2 2 +^gh^h 2 =^gH 2 (A4) 


where subscripts 1 and 2 refer to upstream and down- 
stream conditions, respectively. For given values of Q 
and H, these equations are sufficient to determine u(r), 
h(r), and ly, provided we also specify conditions at r = R. 
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DISCUSSION C. P. Sonett For the distant solar wind the energy density decreases as 1 jr 2 at least. 

At great distance, say a hundred astronomical units, a simple calculation suggests that the 
galactic cosmic ray energy density exceeds that of the solar wind. The cosmic ray gas is, 
of course, an extremely hot gas and it is hard to visualize a supersonic flow in such a gas. 
Now, for the mixture of these two gases (solar wind and cosmic rays) is it still likely that 
one can talk about some kind of shock transition? 

W. I. Ax ford I think it is quite reasonable to treat the two gases separately at least as 
a first approximation, since the coupling is not very strong. One sees plenty of shock 
waves in the interplanetary medium moving through the cosmic ray gas. A strong shock 
does affect the cosmic rays noticeably, but I do not think the shock wave itself is much 
affected. I described these effects briefly at the end of my talk: as the termination of the 
supersonic solar wind flow is approached we expect that the intensity of cosmic rays will 
gradually increase to the interstellar level and then remain roughly constant in the sub- 
sonic region. There may be some acceleration of low energy cosmic rays at the shock as 
discussed by Jokipii [Astrophys. J., 152 , 799 , 1968] , but the more energetic cosmic rays 
do not seem to be noticeably affected. If this acceleration of low energy particles is 
significant and the unmodulated energy density of low energy cosmic rays is much larger 
than observed near the earth (which is possible) then one definitely would have to take 
into account the effect of the cosmic rays on the shock wave. Nevertheless there should 
still be a shock wave of sorts even if it is somewhat smeared out. 

C. P. Sonett It’s true from a structural standpoint, since the gyroradii are larger than 
the solar system, but there are still JXB forces so I’m not sure that one can completely 
ignore the fact that the cosmic ray gas has an enormous temperature. 

Another question is in relation to Geiss’ comment regarding the ratio 4 He/ 3 He. The 
values that he showed for the Apollo samples were 1860-2720, while the value for 
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certain classes of meteorites was about 4000. I wonder, in view of what you said this 
morning, whether the 3 He gradient might be such that the density is higher farther out. 
This class of meteorites have spent most of their lives outside 1 AU and perhaps that would 
account for the concentration. 

W. L Ax ford I think that the 3 He concentration in meteorites could be expected to 
differ from that found in the lunar samples, since the latter are saturated with solar wind 
gases, whereas meteorites have the outer layers (which should also contain solar wind 
material) ablated away. Presumably the 3 He concentration in meteorites is affected by 
cosmic ray implantation and spallation, and any other 3 He deep within the meteorites 
must be primordial in some sense. There is a radial gradient in the density of interstellar 
3 He, but the density within 10 AU is less than that of the solar wind, and of course the 
particle energy is very much less. I doubt that interstellar 3 He could be important as far 
as the interior of meteorites is concerned, and it seems unlikely that it could have much 
effect on the concentrations in lunar dust. 

C. R Sonett Also, there’s the component from charge exchange. 

W. I. Ax ford With the presently accepted values for the density of interstellar gas 
near the sun, I doubt that this could double the 3 He content anywhere in the solar wind 
within 100 AU from the sun. 

J. R. Jokipii I would like to comment on Dr. Sonett’s point and then ask a question. 
It seems to me that if there are regions where the pressure of the cosmic rays is great 
enough, then the cosmic rays will indeed affect the supersonic or the subsonic nature of 
the gas. I’m not sure whether this will ever happen or not. 

W. I. Ax ford Yes, it can in principle (see p.645), but the answer depends upon how 
well the cosmic ray gas is coupled to the solar wind gas, and we have no way of deter- 
mining this without actually exploring the termination region. 

/. R. Jokipii Then one must consider the effect of sound waves actually being 
propagated through the cosmic ray gas. The sound speed is essentially the total pressure 
divided by the density, provided only that the wavelengths are long, enough that the 
cosmic rays and plasma are coupled. We’re talking about 0.01 to 0.1 AU, which is the 
cyclotron radius of a typical cosmic ray and this is small on this scale. Anyway, that’s just 
the comment. The question I have concerns the interstellar density of hydrogen. Presum- 
ably your value is not intended to be an average for the interstellar density of hydrogen. 
It seems to me from some work of Parker and perhaps previous work of others that there 
is trouble with the virial theorem if the density of matter is too low. 

W. I. Ax ford That is a problem and, of course, people have been searching for the 
extra mass for many years. But the mass of neutral hydrogen would not be enough even if 
the number density was as much as 1 cm" 3 . My assumed value of 0.1 cnf 3 is not neces- 
sarily the general density, but it appears to be very close to the local value (within 
~100 pc). As far as the “missing mass” is concerned, there simply is not enough atomic 
hydrogen in the galaxy to matter, and there are observations which suggest that the total 
amount of molecular hydrogen is also not significant in this respect. The only suggestion 
that seems to work Is that the missing mass is in the form of black holes — which has the 
nice feature of being almost impossible to reject on observational grounds. As far as the 
confinement of cosmic rays within the galaxy is concerned, I think we must begin to 
revise our models for the behavior of cosmic rays, perhaps by considering “galactic wind” 
type models in which there is more emphasis on “flow” rather than “leakage” of the cosmic 
ray gas out of the galaxy. 

/. C. Brandt About the mass of the galaxy, that has been missing for decades, I don’t 
think that it’s time to get too worried about it. But there are large areas in the galaxy 
where the number density is apparently 0.1 /cm 3 and I think we just have to live with that. 



M. Dryer I have a question about the term “shock wave.” Whenever we hear the 
term we think of Rankine-Hugoniot conditions, and we happily make comparisons 
between these conditions, including magnetic effects, and observations. Are we still 
talking about the same kind of thin layer in this context, that is, need we worry about 
these things on this scale, or will the meaning of the word “shock wave” be less 
well-specified in this context? 

W. /. Ax ford I don’t see why there should be any problem at all. The scale of the 
system is very large — 100 AU. As long as there is some mechanism which will make the 
shock look “thin” on this scale the Rankine-Hugoniot conditions must be satisfied. Since 
these represent conservation of mass, momentum and energy, there ean hardly be much 
wrong with them. However, one might worry about the effective value of the ratio of 
specific heats, which appears prominently in the Rankine-Hugoniot relations. I doubt that 
the number of degrees of freedom can differ significantly from 3 per particle, but elec- 
tron heat conduction could have a noticeable effect on the transition, as it presumably 
does in the case of shocks observed in the vicinity of the earth. 

C R Sonett Would you care to comment on the time-dependent problem, for 
example, the time-dependent problem that would arise because of the motion of the 
heliosphere through the local arm. This part of the galaxy is very spotty with clouds and 
as the solar system moves through it, the background density varies, the field direction 
changes, and so on. Specifically, how does the information on a change in boundary then 
propagate back into the solar system or heliosphere? 

W. L Axford The “relaxation” time for the system to reach equilibrium is of the 
order of the time the solar wind takes to get to the terminating shock wave. Since the 
solar wind travels an astronomical unit every 4 days, the time scale for the system to 
reach equilibrium is about a year. This is short compared with any possible time scale for 
changes occurring outside the solar system. At a number density of 0.1 cm" 3 , the mean 
free path is about IQ 17 cm. At 20 km/ sec, which is roughly the velocity of the interstellar 
gas past the sun, an interval of 5.10 10 sec is required for a “cloud” with a scale size equal 
to one mean free path to go past the solar system. So a cloud which is only one mean free 
path across takes about 2000 years to go past the sun at this speed. If you take 
1000 mean free paths for the size of an interstellar cloud, one sees that the heliosphere 
should be stationary over time scales of a million years. 

C. j P. Sonett Is that true also for the wave field? Supposing there were a wave field 
on the outside, as we suspect there probably is. 

W. I Axford Well, the waves that I would imagine to be important are Alfven waves. 
The Alfven speed is also quite small, something like 30 km/ sec, so the time scale is still 
going to be of the order of 10 6 years if the wavelengths are comparable to cloud 
dimensions. 



